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GLOBAL  MAPS  OF  foF2  DERIVED  FRuM 
OBSERVATIONS  AND  THEORETICAL  VALUES 


Charles  M.  Rush* 
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Frank  G.  Stewart* 

Rita  Reasoner* 

Observations  of  the  F2-region  critical  frequency,  foF2,  and 
values  determined  from  the  time -dependent  continuity  equation  for 
ions  and  electrons  in  the  ionosphere  have  been  used  to  develop  a 
new  set  of  numerical  coefficients  to  represent  the  global  vari¬ 
ation  of  foF2.  The  coefficients  developed  in  this  study,  like 
those  in  earlier  investigations,  permit  monthly  median  hourly 
values  of  foF2  to  be  obtained  at  any  location  around  the  globe  for 
any  month  of  the  year  and  solar  activity  level.  The  procedures 
undertaken  to  develop  the  new  numerical  coefficients  are  described 
in  this  report.  Comparisons  between  foF2  determined  using  older 
sets  of  numerical  coefficients  and  foF2  determined  using  the  new 
set  of  coefficients  are  given  along  with  a  description  of  how  well 
each  set  of  coefficients  specifies  and  predicts  the  observed  vari¬ 
ations  in  the  F2-region  critical  frequency. 

Key  Words:  ionospheric  maps;  critical  fre^ency;  time-dependent 
continuity  equation 


1.  INTRODUCTION 

The  F2  region  of  the  ionosphere  exhibits  a  profound  influence  upon  the 
performance  of  long  distance  HF  propagation  circuits  and  upon  radio  signals 
(HF  and  VHF)  that  propagate  through  the  ionosphere.  The  structure  of  the  F2 
region  varies  with  time  and  space  and  imparts  comparable  variations  upon  the 
radio  signals  that  traverse  it.  The  critical  frequency  of  the  F2  region. 
foF2,  is  a  fundamental  parameter  that  is  used  in  the  specification  and  predic¬ 
tion  of  the  structure  of  the  F2  region.  Global  maps  of  foF2  are  used  as  the 
oasis  for  various  ionospheric  models  (Nisbet,  1971),  HF  propagation  prediction 
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methods  (CCIR,  1982a;  Teters  et  al.,  1983),  and  transionospheric  propagation 
models  (Bent  et  al.,  1978).  These  maps  (denoted  hereafter  as  CCIR  maps)  are 
generated  from  numerical  coefficients  obtained  by  performing  a  spherical  har¬ 
monic  analysis  on  observed  monthly  median  values  of  foF2  described  by  Jones 
and  Gal  let  (1962).  The  accuracy  of  the  global  representation  of  foF2  obtained 
from  the  CCIR  maps  is  dependent  upon  the  geographical  distribution  of  the 
stations  whose  data  were  used  in  the  generation  of  the  numerical  coeffi¬ 
cients.  Where  observations  were  available  for  inclusion  in  the  analysis  that 
generated  the  coefficients,  the  accuracy  of  foF2  is  quite  reasonable  (King  and 
Slater,  1973).  At  locations  where  data  were  not  available  (such  as  oceans), 
the  accuracy  of  foF2  determined  from  the  coefficients  is  questionable. 

In  a  recent  paper.  Rush  et  al.  (1983)  described  how  theoretically  derived 
values  of  the  F2-region  critical  frequency,  foF2,  could  be  used  to  improve  the 
global  representation  of  foF2  that  is  achievable  using  the  CCIR  numerical 
coefficients  (CCIR,  1982b).  The  improvement  is  obtained  by  combining  values 
of  foF2  derived  from  the  time-dependent  continuity  equation  for  ions  (and 
electrons)  with  monthly  median  observations  of  foF2  to  obtain  data  bases  that 
can  be  subjected  to  spherical  harmonic  analysis.  The  harmonic  analysis  then 
yields  coefficients  that  represent  the  global  variation  of  foF2  that  is  con¬ 
sistent  with  both  observations  and  the  current  understanding  of  the  physical 
mechanisms  responsible  for  the  large-scale  ionospheric  variations. 

The  results  presented  by  Rush  et  al.  (1983)  were  confined  to  the  months 
of  July,  September,  and  December  and  the  years  1975  and  1978.  The  results  for 
the  year  1975  were  taken  as  being  indicative  of  sunspot  minimum  conditions  and 
those  for  1978  were  indicative  of  sunspot  maximum  conditions.  A  limited  num¬ 
ber  of  comparisons  between  foF2  values  derived  from  the  CCIR  maps  and  the  foF2 
values  derived  from  the  new  coefficients  were  also  given.  Since  that  paper 
was  published,  numerical  coefficients  for  each  month  of  a  sunspot  minimum  and 
a  sunspot  maximum  year  have  been  developed  that  are  based  upon  the  same  prin¬ 
ciples  cs  given  by  Rush  et  al.  (1983).  It  is  the  purpose  of  this  report  to 
describe  the  most  recent  results  obtained  and  to  present  detailed  comparisons 
of  observed  foF2  values  with  values  of  foF2  derived  from  the  CCIR  maps  and 
with  values  of  foF2  derived  from  the  new  coefficients.  For  the  sake  of 
completeness,  the  results  obtained  for  the  months  of  July,  September,  and 
December  1975  and  1978  which  were  discussed  by  Rush  et  al .  (1983)  are  also 
presented  here. 
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In  the  next  section,  the  data  and  the  sources  of  data  used  in  the  devel¬ 
opment  of  the  new  foF2  coefficients  are  described.  Observations  ot  monthly 
median  hourly  values  of  foF2  obtained  from  vertical  incidence  ionosondes 
located  throughout  the  world  were  used  as  the  primary  data  source.  In  the 
third  section,  the  use  of  the  time-dependent  continuity  equation  to  derive 
values  of  foF2  is  presented.  In  Section  4,  the  analysis  procedures  employed 
to  produce  the  new  coefficients  are  briefly  described;  and,  in  Section  5,  the 
new  maps  are  described  and  the  results  of  the  many  comparisons  that  were 
undertaken  as  part  of  this  study  are  presented.  In  the  last  section,  the 
overall  results  are  summarized  and  the  conclusions  drawn  from  them  are  given. 


2.  DATA  AND  DATA  SOURCES 

In  this  study,  numerical  coefficients  that  represent  the  global  varia¬ 
tions  of  foF2  have  been  determined  for  each  month  for  July  1975  to  June  1976 
and  for  July  1978  to  June  1979.  The  coefficients  for  July  1975  to  June  1976 
are  representative  of  solar  minimum  conditions  and  those  for  July  1978  to  June 
1979  are  representative  •  :  olar  maximum  conditions.  Table  1  lists  the  12- 
month  average  Zurich  sunspot  number  for  each  of  the  months  for  which  coeffi¬ 
cients  were  derived.  It  is  evident  from  the  table  that  for  solar  minimum 
(1975/1976)  the  sunspot  numbers  do  not  vary  greatly  from  month-to-month  and 
are  typically  on  the  order  of  10  to  15  units.  For  solar  maximum  (1978/1979), 
however,  the  sunspot  number  increases  continuously  from  95  in  July  1978  to  154 
in  June  1979.  Thus,  the  coefficients  obtained  in  this  study  for  solar  maximum 
may  not  accurately  represent  the  highest  solar  maximum  conditions  that  are  ob¬ 
served,  particularly  for  the  months  July  through  December  1978. 

In  developing  the  first  set  of  numerical  maps  of  ionospheric  data,  Jones 
and  Gal  let  (1962)  found  it  necessary  to  ensure  that  foF2  data  existed  at  suf¬ 
ficient  locations  so  that  the  mapped  fields  did  not  have  negative  values. 

These  negative  values  resulted  from  instabilities  in  the  maps  in  areas  where 
data  were  either  nonexistent  or  sparse.  In  order  to  prevent  instabilities, 
Jones  and  Gal  let  employed  a  "screen  analysis"  whereby  observations  available 
at  specific  locations  were  extrapolated  to  areas  devoid  of  data  taking  account 
of  the  geomagnetic  field  control  of  the  F2  region.  Inherent  in  this  procedure 
is  the  assumption  that  the  temporal  variations  in  foF2  at  one  location  (obser¬ 
vation  point)  are  simply  related  to  those  at  another  location  (extrapolated 
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Table  1.  Smoothed  Observed  Sunspot  Numbers  for 
July  1975  Through  June  1976  and 
July  1978  Through  June  1979 


i 


JAN  FEB  MAR  APR  MAY  JUN 

JUL  AUG  SEP  OCT  NOV  DEC 

i  1975 

15.0  14.3  14.5  15.6  16.3  16.5 

|1976  15.2  13.2  12.2  12.6  12.5  12.2 

1 

| 

S 1978 

95.0  104.0  108.4  111.0  113.3  116.7 

1 1979  122.8  130.4  136.1  141.0  147.4  153.7 

point)  in  a  relatively  simple  manner.  This  assumption  may  be  valid  for  regions 
over  which  extrapolation  does  not  have,  to  be  extended  to  great  distances  (see 
for  example  Rush,  1976).  However,  the  inaccuracies  in  the  maps  over  the  ocean 
areas,  for  example,  point  up  the  need  for  a  better  estimate  of  the  behavior  of 
foF2  in  regions  of  the  earth  inaccessible  to  routine  observations. 

In  the  development  of  the  new  set  of  numerical  coefficients  discussed  in 
this  report,  the  screen  analysis  of  Jones  and  Gal  let  (1962)  was  abandoned  in 
favor  of  a  more  physically  based  approach.  In  order  to  obtain  a  uniform  data 
set  that  could  be  subjected  to  the  spherical  harmonic  analysis  techniques  used 
by  Jones  and  Gallet,  observations  of  foF2  obtained  from  vertical -incidence 
ionosondes  were  combined  with  foF2  values  determined  from  the  time-dependent 
continuity  equation  for  ions  (and  electrons).  In  this  section,  the  values  of 
foF2  determined  directly  from  observations  or  inferred  from  observations  will 
be  discussed.  The  values  of  foF2  determined  by  use  of  the  continuity  equa¬ 
tions  will  be  discussed  in  the  next  section. 

Monthly  median  observations  of  r"oF2  obtained  from  the  stations  given  in 
Table  2  form  the  basis  for  the  spherical  harmonic  analysis.  In  the  table,  the 
station  name  and  location  in  geographic  coordinates  are  given.  Also  given  in 
the  table  is  the  letter  A  or  B  under  each  of  the  two  solar  cycle  epochs.  The 
letter  A  refers  to  data  that  were  observed  at  the  station  for  the  time  period 
1975/1976  and/or  1978/1979.  The  letter  B  refers  to  data  observed  at  tne  sta¬ 
tion  in  question  for  a  different  solar  epoch  and  adjusted  to  the  solar  epochs 
in  1975/1976  or  1978/1979  using  a  polynomial  interpolation  procedure 


Table  2.  lonosonde  Stations  Providing  Data  Used 

in  the  Determination  of  foF2  Coefficients 


Station  Name 

Latitude 

(°N) 

Longitude 

(°E) 

1975/1976 

1978/1979 

Arctic  I  (USSR)  Floating 

Ice 

Island 

B 

B 

Arctic.  II  (USSR)  Floating 

Ice 

Island 

B 

B 

Arctic  III  (USSR)  Floating 

Ice 

Island 

B 

B 

Heiss  Island,  USSR 

80.6 

58.0 

A 

B 

Thule,  Greenland 

76.4 

291.3 

B 

B 

Resolute  Bey,  Canada 

74.7 

265.1 

A 

A 

Dixon  Island,  USSR 

73.5 

80.4 

B 

B 

Barrow,  Alaska,  USA 

71.3 

203.2 

B 

B 

Godhavn,  Greenland 

69.3 

306.5 

B 

B 

Kiruna,  Sweden 

67.8 

20.4 

A 

A 

Sodanxyla,  Finland 

67.4 

26.6 

A 

A 

Salekhard,  USSR 

66.5 

66.7 

A 

A 

Lycksele,  Sweden 

64.7 

18.8 

A 

A 

Archangel sk,  USSR 

64.6 

40,5 

- 

A 

Providenya,  USSR 

64.4 

186.6 

B 

B 

Tunguska,  USSR 

61.6 

90.0 

- 

A 

Yakutsk,  USSR 

62.0 

129.6 

- 

A 

Narssarscuaq,  Greenland 

61.2 

314.6 

B 

B 

Anchorage,  Alaska,  USA 

61. 2 

210.1 

B 

B 

Leningrad,  USSR 

60.0 

30.7 

A 

A 

Nurmijarvi,  Finland 

60.5 

24.6 

A 

A 

Uppsala,  Sweden 

59.8 

17.6 

A 

A 

Churchill.  Canada 

58.8 

265.8 

A 

A 

Gorky,  USSR 

56.1 

44.3 

A 

A 

Sverdlovsk,  USSR 

56.7 

61.1 

A 

A 

Tomsk,  USSR 

56.5 

84.9 

A 

A 

Uul iusruh/Rugen,  Federal 

54.6 

13.4 

A 

A 

Republic  of  Germany 

Moscow,  USSR 

55.5 

37.3 

A 

A 

Goosebay,  Canada 

53.3 

299.2 

A 

A 

De  Bilt,  Netherlands 

52.1 

5.2 

A 

A 

Pruhonice,  Czechoslovakia 

50.0 

14.6 

A 

A 

Miedzeszyn,  Poland 

52.2 

21.2 

A 

A 

Irkutsk,  USSR 

52.5 

104.0 

A 

A 

Slough,  England 

51.5 

359.4 

A 

A 

Adak,  Alaska,  USA 

51.9 

183.  * 

B 

B 

Lindau,  Federal  Republic  of 

51.6 

10.1 

A 

A 

Ge rmany 

Dourbes,  Belgium 

50.1 

4.6 

A 

A 

Rostov,  USSR 

47.2 

39.7 

A 

A 

Winnipeg,  Canada 

49.8 

265.6 

A 

Freiburg,  Federal  Republic 

of 

48.1 

7.6 

A 

- 

Germany 

Khabarovsk,  USSR 

48.5 

125.1 

A 

A 

Lannion,  France 

48.5 

356.7 

A 

A 

Budapest,  Hungary 

46.7 

21.2 

\ 

<• 

Saint  John's,  Canada 

47.6 

307.3 

A 

- 

Table  2.  Ionosonde  Stations  Providing  Data  Used  in  the 
Determination  of  foF2  Coefficients  (Continued) 


Station  Name 

Latitude 

(°N) 

Longitude 

(°E) 

1975/1976 

1978/1979 

[• 

Poitiers,  France 

46.6 

0.3 

A 

A 

Wakkanai,  Japan 

45.4 

141.7 

A 

A 

Ottawa,  Canada 

45.4 

284.1 

A 

A 

Alma-Ata,  USSR 

43.2 

76.9 

A 

A 

.  ■  . .  ,  *  •  * 

Tbilisi,  USSR 

41.7 

44.8 

A 

A 

r* 

Rome,  Italy 

41.8 

12.5 

A 

B 

/’•  > 

Tashkent,  USSR 

41.3 

69.6 

- 

A 

'  •  Jl,' 

Tortosa,  Spain 

40.8 

0.3 

A 

B 

Boulder,  Colorado,  USA 

40.0 

254.7 

A 

A 

Akita,  Japan 

39.7 

140.1 

A 

A 

—•  ■  — 

Gibilmanna,  Sicily 

38.0 

14.0 

- 

A 

Ashkabad,  USSR 

37.9 

58.3 

A 

A 

, 

Seoul ,  South  Korea 

37.4 

127.0 

A 

A 

Wallops,  Virginia,  USA 

37.9 

284.5 

A 

A 

Point  Arguello,  California,  USA 

35.6 

239.4 

A 

A 

-.’V 

Tokyo,  Japan 

35.7 

139.5 

A 

A 

White  Sands,  New  Mexico,  USA 

32.3 

253.5 

- 

A 

Yamagawa,  Japan 

31.2 

130.6 

A 

A 

•* 

Wuhan,  China 

30.6 

114.4 

A 

A 

Delhi,  India 

28.6 

77.2 

B 

B 

Okinawa,  Japan 

26.3 

127.8 

A 

A 

Taipei ,  Taiwan 

25.0 

121.2 

A 

A 

Ahmedabad,  India 

23.0 

72.6 

B 

B 

;  • 

Hong  Kong,  Southeast  Asia 

22.3 

114.2 

B 

E 

.  ■ ",  • 1 ,  •* 

Maui,  Hawaii,  USA 

20.8 

203.5 

A 

A 

•  v 

Mexico  City,  Mexico 

19.4 

260.3 

A 

A 

p*  /• . 

Dakar,  Senegal 

14.7 

342.6 

A 

B 

Manila,  Phillipines  Island 

14.7 

121.1 

A 

A 

Bangkok,  Thailand 

13.7 

100.6 

A 

A 

r.ii .. 

Ouagadougou,  Upper  Volta 

12.4 

358.5 

A 

A 

*  Vt~,  . 

Djibouti,  Northeast  Africa 

11.6 

42.8 

A 

A 

•• '«  . 

Kodaikanal,  India 

10.2 

77.5 

A 

- 

Thumba,  India 

8.6 

76.9 

A 

- 

Vanimo,  New  Guinea 

-2.7 

141.3 

A 

A 

Huancayo,  Peru 

-12.0 

284.7 

A 

A 

• 

Tahiti ,  Paci fic  Ocean 

-17.7 

210.7 

A 

A 

■  • 

Tsumeb,  Namibia 

-19.2 

17.7 

A 

- 

Townsvil le,  Austral  ia 

-19.3 

146.7 

A 

A 

Raratonga,  Cook  Island 

-21.2 

200.2 

A 

A 

**• 

Johannesburg,  South  Africa 

-26.1 

28.1 

A 

A 

v  v 

Tucuman,  Argentina 

-26.9 

294.6 

A 

A 

*  ' 

Brisbane,  Australia 

-27.5 

152.9 

A 

A 

Norfolk  Island,  Pacific  Ocean 

-29.0 

168.0 

A 

A 

.  •  ,  • . * 

Mundaring,  Australia 

-32.0 

116.2 

A 

A 

*  *  ' 

Capetown,  South  Africa 

-34.1 

18.3 

A 

- 

'■  y%. 

Buenos  Aires,  Argentina 

-34.5 

301.5 

A 

- 

Canberra,  Australia 

-35.3 

149.0 

A 

A 

• 
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Table  2.  Ionosonde  Stations  Providing  Data  Used  in  the 
Determination  of  foF2  Coefficients  (Concluded) 


Station  Name 

Latitude 

r  n) 

Longitude 

.{'El. 

1975/1976 

1978/1979 

Concepcion,  Chile 

-36.6 

287.0 

A 

A 

Auckland,  New  Zealand 

-37.0 

175.0 

A 

A 

Hobart,  Australia 

-42.9 

147.2 

A 

A 

Christchurch,  New  Zealand 

-43.6 

172.8 

A 

A 

Kerguelen,  Indian  Ocean 

-49.4 

70.3 

A 

A 

Port  Stanley,  Falkland  Islands 

-51.7 

302.2 

A 

A 

Campbell  Island,  Pacific  Ocean 

-52.5 

169.2 

A 

A 

South  Georgia,  Atlantic  Ocean 

-54.3 

323.5 

A 

- 

Argentine  Island,  Antarctica 

-65.2 

295.7 

A 

A 

Casey  Base,  Antarctica 

-66.2 

110.5 

B 

B 

Mirny,  Antarctica 

-66.5 

93.0 

B 

B 

Terre  Adelie,  Antarctica 

-66.7 

140.0 

A 

A 

Mawson,  Antarctica 

-67.6 

62.9 

B 

A 

Syowa  Base,  Antarctica 

-69.0 

39.6 

A 

A 

Sanae,  Antarctica 

-70.3 

357.6 

A 

A 

Halley  Bay,  Antarctica 

-75.5 

333.4 

A 

A 

Scott  Base,  Antarctica 

-77.8 

166.8 

A 

A 

Vostok,  Antarctica 

-78.4 

106.9 

B 

B 

Byrd  Station,  Antarctica 

-80.0 

240.0 

B 

B 

South  Pole,  Antarctica 

-90.0 

0.0 

B 

B 

Key:  The  symbol  A  indicates  that  median  observations  for  the  time  period  in 
question  were  used  in  the  analysis. 


The  symbol  B  indicates  that  observations  were  used  to  predict  the  median 
value  of  foF2  for  the  time  period  in  question. 


{Crow  and  Zacharisen,  1960).  To  obtain  values  of  foF2  from  the  Interpolation 
procedure  for  a  given  location,  observations  at  that  location  had  to  be  avail¬ 
able  for  the  month  in  question  for  other  solar  epochs.  Data  of  this  type  are 
referred  to  in  this  report  as  8  data.  The  use  of  data  similar  to  B  data  in 
other  geophysical  disciplines  has  been  pointed  out  by  Rush  et  al.  (1982). 

Actual  observed  monthly  median  data  are  denoted  as  A  data. 

Figure  1  depicts  the  distribution  of  A-  and  B-data  sources  throughout  the 
world.  The  figure  vividly  illustrates  the  lack  of  observed  data  spread  uni¬ 
formly  throughout  the  globe.  Obviously  A-  and  B-data  sources  are  available  over 
land-based  regions.  Most  of  the  globe,  however,  is  not  covered  with  these 
observations.  This  situation  was  ameliorated  by  using  the  time-dependent 
continuity  equation  to  generate  foF2  values  in  the  data-sparse  areas.  In  actual 
fact,  the  continuity  equation  was  used  to  determine  foF2  values  at  middle  and 
high  latitudes  (latitudes  greater  than  about  25°  north  or  south  magnetic 
latitude). 

The  data  needed  for  the  unobserved  regions  in  the  equatorial  ionosphere 
were  not  derived  in  as  rigorous  a  fashion  as  those  for  the  unobserved  regions  at 
other  latitudes.  This  is  due  primarily  to  the  fact  that  for  the  time-deoendent 
continuity  equation  to  yield  realistic  values  of  foF2  at  low  magnetic  latitudes, 
the  vertical  electrodynamic  drift  at  the  equator  must  be  included  in  the  calcu¬ 
lation.  This  results  in  calculations  that  are  very  lengthy  and  time-consuming 
even  for  modern  computers.  It  was  determined  that  it  was  impractical  in  this 
investigation  to  obtain  values  of  foF2  needed  for  the  low  latitudes  from  the 
theoretical  or  time-dependent  continuity  equation.  Rather,  what  was  done  was 
to  take  foF2  observations  at  the  available  low-latitude  locations  and  plot 
them  onto  global  maps  for  each  hour  of  local  time.  The  variation  of  foF2 
plotted  at  constant  local  time  is  much  smoother  than  when  plotted  in  constant 
universal  time.  Contours  of  constant  foF2  were  then  constructed  at  each  hour 
of  local  time  for  each  month  from  the  plotted  values  of  foF2.  It  was  found 
that  distinct  and  noticeable  differences  in  the  shape  and  form  of  the  contours 
had  to  be  drawn  for  regions  of  the  globe  where  the  magnetic  field  has  an  east¬ 
ward  declination  and  where  it  has  a  westward  declination.  Once  the  contours 
were  drawn  for  each  hour  of  local  time  for  a  given  month,  the  local  time  maps 
were  used  to  determine  foF2  at  specific  points  in  the  low-latitude  ionosphere. 
These  low-'^itude  values  were  then  combined  with  the  values  of  foF2  deter¬ 
mined  fror  the  time-dependent  continuity  equations  and  the  A  and  8  data  to 
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provide  the  data  base  needed  to  generate  numerical  coefficients  from  the 
spherical  harmonic  analysis  of  Jones  and  Gal  let  (1962). 

3.  THEORETICALLY  DERIVED  DATA 

As  mentioned  in  previous  sections,  the  time-dependent  ion  continuity 
equation  was  used  to  calculate  values  of  foF2  at  middle  and  high  latitude 
regions  where  observations  of  foF2  were  not  available.  The  values  of  foF2 
thus  determined  are  referred  to  hereafter  as  theoretical  values.  They  were 
obtained  by  use  of  the  time-dependent  ion  continuity  equation  as  described  by 
Anderson  (1973)  and  applied  to  ionospheric  mapping  studies  by  Rush  et  al . 
(1983). 

The  continuity  equation  is  given  by 

3N. 

-jj-.V.  (N»t)  -L,  (1) 

where  Ni  is  the  ion  density;  Pj,  the  ion  production  rate,  Li,  the  loss  rate; 
and  Vi,  the  transport  velocity.  In  the  ionosphere,  plasma  is  transported 
along  the  geomagnetic  field  lines  by  diffusion  and  neutral  winds  and  perpen¬ 
dicular  to  the  field  lines  primarily  by  electrodynamic  (E  x  B)  drifts.  For 
the  purposes  of  this  investigation,  it  was  assumed  at  middle  and  high  lati¬ 
tudes  that  the  transport  of  plasma  along  the  field  lines  greatly  exceeded  the 
transport  across  the  field  lines.  The  movement  of  plasma  due  to  electro¬ 
dynamic  drift  was,  therefore,  neglected  permitting  a  much  simpler  and  computa¬ 
tionally  efficient  solution  to  Equation  (1). 

The  continuity  equation  was  solved  numerically  following  Anderson  (1973), 
to  give  (=  Nfi,  electron  density)  as  a  function  of  altitude,  latitude,  and 
local  time.  In  arriving  at  a  solution  for  it  was  assumed  that  the  only 
ion  of  concern  was  atomic  oxygen  (0+),  a  valid  approximation  for  F2-region 
heights.  The  parameters  needed  to  solve  the  ion  continuity  equation  were 
obtained  from  models  of  the  neutral  composition  and  temperature,  ion  and 
electron  temperatures,  production,  loss  and  diffusion  rates,  and  neutral  wind 
models.  These  models  were  determined  for  the  appropriate  month  and  sunspot 
number  given  in  Table  1.  In  addition,  a  model  of  the  geomagnetic  field  had  to 
be  assumed  in  order  to  arrive  at  a  solution  to  Equation  (I).  The  details  of 
the  models  of  the  parameters  used  to  develop  theoretical  values  of  foF2  are 
discussed  below. 
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where 


0O2  =  3.3  x  10*®  cm"*  sec"*, 
oNg  =  3.9  x  10*®  cm"*  sec"*, 
and  oO  =  3.7  x  10*®  cm'*  sec"*. 

The  terms  [O2],  [N2],  and  [0]  represent  the  number  densities  of  molecular 
oxygen,  molecular  nitrogen,  and  atomic  oxygen,  respectively.  These  values, 
along  with  Tn,  were  derived  from  the  neutral  atmosphere  model  and  vary  with 
time  in  accord  with  the  MSI S  model. 


Table  3.  Rate  of  Ion  Production,  P(sec"*) 

1975/1976 _ 1978/1979 


[ 

!  Month 

Northern 

Hemishpere 

Southern 

Hemispshere 

Northern 

Hemisphere 

Southern  ] 
Hemisphere  !• 

[  July 

2.3 

X 

10"7 

2.4 

X 

10-7 

4.0 

X 

10"7 

4.0 

X 

10"7  j 

-  August 

2.3 

X 

10"7 

2.4 

X 

10’7 

4.0 

X 

10"7 

4.0 

X 

10"7  i 

■  September 

2.3 

X 

10'7 

2.3 

X 

10"7 

4.0 

X 

lO"7 

4.0 

X 

10"7  1 

October 

2.3 

X 

10~7 

2.3 

X 

10"7 

4.5 

X 

10'7 

4.5 

X 

10"7  ; 

November 

2.3 

X 

10~7 

2.3 

X 

10"7 

4.5 

X 

10"7 

4.5 

X 

10" 7  : 

December 

2.3 

X 

10"7 

2.3 

X 

10"7 

4.5 

X 

10'7 

4.5 

X 

10"7  . 

January 

2.3 

X 

10"7 

2.3 

X 

I0"7 

5.0 

X 

K-* 

0 

t 

4.5 

X 

10"7 

February 

2.3 

X 

10'7 

2.3 

X 

iir7 

5.0 

X 

10"7 

4.5 

X 

10"7  ' 

March 

2.5 

X 

10"7 

2.5 

X 

10“7 

6.0 

X 

10"7 

6.0 

X 

10"7  • 

April 

2.3 

X 

1G'7 

2.3 

X 

I0"7 

6.0 

X 

10"7 

6.0 

X 

10"7  : 

May 

2.3 

X 

I0"7 

2.4 

X 

10-7 

5.0 

X 

10-7 

6.0 

X 

10-7  . 

June 

2.3 

X 

10"7 

2.3 

X 

iq-? 

4.0 

X 

10-7 

4.0 

X 

10"7 

_ 1 

3.4  Neutral  Wind  Models 

Values  of  meridional  and  zonal  neutral -air  winds  for  each  of  the  months 
were  deduced  by  fitting  theoretically  calculated  values  of  foF2  to  actual 
monthly  median  observations  at  specific  aid -latitude  locations.  The  conti¬ 
nuity  equation  was  used  to  calculate  hourly  values  of  foF2  for  a  location 
where  foF2  observations  were  available  and  the  meridional  component  of  the 
wind  was  adjusted  so  that  agreement  between  observed  and  calculated  values 
were  obtained.  The  calculation  was  then  repeated  for  a  location  where  the 
foF2  data  were  available  and  where  the  magnetic  latitude  was  essentially  the 
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same  as  for  the  first  calculation  but  the  declination  of  the  geomagnetic  field 
was  different.  Using  the  same  meridional  wind,  the  zonal  wind  component  was 
adjusted  so  that  calculated  values  of  foF2  agreed  with  observations  at  the 
second  location. 

Meridional  and  zonal  neutral  wind  models  were  determined  for  each  of  the 
months  in  the  1975/1976  and  1978/1979  time  periods.  Different  models  for  the 
same  months  were  obtained  for  the  northern  and  southern  hemispheres  of  the 
globe,  and  the  winds  were  assumed  to  be  independent  of  altitude. 

Figure  2  shows  the  neutral  wind  models  for  the  months  of  January  1976  and 
1979  for  both  the  northern  and  southern  hemispheres.  The  variations  seen  in 
the  zonal  and  meridional  winds  in  Figure  2  are  typical  of  the  neutral  wind 
patterns  observed  for  all  the  months  used  in  this  analysis.  The  zonal  wind  is 
directed  eastward  during  the  evening  hours  and,  for  some  months,  in  the  hours 
immediately  following  midnight.  During  solar  maximum,  the  eastward  zonal  wind 
prevails  in  the  winter  hemisphere  around  the  sunrise  hours.  This  is  not  the 
case  at  solar  minimum  when  the  zonal  wind  is  either  absent  or  directed  to  the 
west  during  the  sunrise  hours.  The  meridional  wind  is  directed  equatorward 
during  the  night  ar.d  early  morning  hours  and  is  either  absent  or  directed 
toward  the  poles  during  the  daylight  hours. 

The  wind  patterns  for  each  month  used  in  this  investigation  are  shown  in 
figures  in  Appendix  A.  Certain  generalities  can  be  made  about  the  variations 
in  the  wind  patterns  by  studying  the  figures  in  Appendix  A.  The  eastward 
zonal  wind  is  larger  in  the  nighttime  winter  hemisphere  than  in  the  summer 
hemisphere  during  solar  minimum,  but  not  so  during  solar  maximum.  Except  for 
the  equinoctial  months,  the  eastward  directed  zonal  wind  during  solar  minimum 
is  larger  than  during  solar  maximum.  The  meridional  winds  display  similar 
behavior.  The  equatorward  meridional  wind  in  the  nighttime  winter  hemisphere 
is  generally  greater  than  that  in  the  summer  hemisphere.  The  daytime  poleward 
meridional  wind  also  tends  to  be  larger  during  the  winter  months  than  during 
the  sumner  months. 

The  neutral -air  winds  at  F2-region  altitudes  result  from  the  temperature 
and  pressure  gradients  set  up  in  the  neutral  atmosphere  by  diurnal  heating  and 
cooling  (King  and  Kohl,  1965).  The  speed  of  the  winds  is  slowed  by  colli- 
sional  interactions  between  the  neutral  atmosphere  and  the  ions  at  the 
F2-region  heights.  During  times  when  the  ionization  density  is  large  (such  as 
solar  maximum  compared  to  solar  minimum),  the  collisior.al  forces  or  "ion-drag" 
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forces  reduce  the  speed  of  the  wind.  The  general  variations  In  the  wind  pat¬ 
terns  and  the  magnitude  of  the  winds  shown  in  Appendix  A  are  consistent  with 
the  MSIS  neutral  atmospheric  model  and  the  diurnal,  seasonal,  and  solar  cycle 
variations  of  the  ionization  density  in  the  F2  region. 

In  the  F  region,  the  neutral  wind  moves  ionization  along  geomagnetic 
field  lines.  The  equatorward  meridional  wind  moves  ionization  up  in  altitude, 
and  the  poleward  meridional  wind  moves  ionization  downward.  The  zonal  wind 
imparts  an  upward  or  downward  drift  to  the  ionization  depending  upon  the  rela¬ 
tive  direction  between  the  wind  velocity  and  the  declination  of  the  geomag¬ 
netic  field.  In  the  southern  hemisphere  where  the  inclination  of  the  geomag¬ 
netic  field  is  directed  upward,  an  eastward  wind  will  impart  a  downward  com¬ 
ponent  to  the  ionization  when  the  field  is  declined  westward  and  an  upward 
component  when  it  is  declined  eastward.  A  westward  wind  in  the  southern 
hemisphere  gives  rise  to  an  upward  component  when  declination  is  westward  and 
a  downward  component  when  the  declination  is  eastward.  The  effect  in  the 
northern  hemisphere  is  just  the  opposite  because  the  inclination  of  the 
magnetic  field  lines  is  directed  downward.  Figure  3  provides  a  simplified 
illustration  of  the  movement  of  ionization  along  the  geomagnetic  field  due  to 
zonal  winds  in  the  southern  hemisphere  and  northern  hemisphere  for  east  and 
west  geomagnetic  declinations. 

The  interaction  between  the  neutral -air  wind  zonal  component  and  the  geo¬ 
magnetic  field  will  give  rise  to  very  different  diurnal  variations  of  foF2  at 
eastward  and  westward  declinations.  Simply  extrapolating  observation*  from 
one  location  to  another  and  taking  into  account  only  the  magnetic  latitude 
dependence  of  foF2  could  result  in  placing  data  in  areas  where  the  declination 
is  vastly  different.  This  would  give  rise  to  foF2  values  in  inaccessible 
regions  that  are  unrealistic  and  unrepresentative  of  the  area.  The  magnitude 
of  the  error  made  will  depend  upon  the  difference  between  the  declination  at 
the  observation  location  and  the  extrapolated  location.  Figure  4,  which  shows 
the  magnetic  declination  of  the  1975  epoch  determined  by  the  United  States 
Defense  Mapping  Agency,  illustrates  the  problem  of  simple  extrapolation. 

Moving  observations  from  land  areas  anywhere  in  the  world  to  adjacent  ocean 
areas  results  in  placing  the  observations  in  regions  with  vastly  different 
magnetic  declinations.  By  using  the  time -dependent  continuity  equation  to 
derive  values  of  foF2  at  inaccessible  regions,  the  magnetic  inclination  and 
declination  dependencies  in  the  F2  region  are  properly  taken  into  account. 
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geomagnetic  field  epoch. 


3.5  Geomagnetic  Field  Model 

A  realistic  geomagnetic  field  model  was  used  to  transform  the  spherical 
polar  coordinate  system  r,  6,  and  <j>  to  one  describing  directions  parallel  and 
perpendicular  to  B.  The  "real"  field  is  given  by  Euler  potentials  a  and  $ 
where  B  =  Va  x  VB  (Stern,  1965,  1967),  and  a  full  description  of  the  technique 
is  covered  in  the  paper  by  Anderson  (1973). 

3.6  Solution  to  the  Full  Continuity  Equation 

The  time-dependent  continuity  equation  was  solved  numerically  using  a 
Crank-Nicolson  (Crank  and  Nicolson,  1947)  implicit  finite  differencing  scheme. 
The  two  boundaries  are  at  125  km  altitude  at  the  northern  and  southern  ends  of 
the  field  line.  Boundary  conditions  for  the  ion  density  are  photochemical 
equilibrium  during  the  day  and  10  ions/cm  at  night.  However,  the  boundary 
conditions  are  not  critical,  as  the  appropriate  ion  densities  are  reached 
within  a  few  space  steps.  The  changing  size  of  the  magnetic  flux  tube  is 
taken  into  account  in, the  calculations. 

The  importance  of  using  values  of  foF2  that  are  derived  from  the  conti¬ 
nuity  equation  to  supplement  observations  needed  to  develop  realistic  iono¬ 
spheric  maps  can  be  seen  in  Figure  5.  The  figure  shows  the  diurnal  variations 
of  foF2  for  September  1978  at  three  different  longitudes  where  the  declination 
of  the  magnetic  field  line  changes  but  the  dip  angle  is  approximately  the 
same.  [This  figure  is  the  same  as  Figure  4  ir.  Rush  et  al .  (1982).]  During 
daytime  hours,  the  westward  declination  field  line  has  the  greatest  peak 
electron  density  because,  in  the  southern  hemisphere,  the  early  morning, 
westward-blowing  zonal  wind  (see  Figure  A. 12)  at  50  m/sec  transports  ioni¬ 
zation  upward  from  below  where  the  production  rate  is  greater  to  higher  alti¬ 
tudes  v-here  the  loss  rate  is  less.  In  geographic  regions  where  the 
declination  is  nearly  zero,  the  effect  of  the  zonal  wind  is  minimal  and  the 
results  shown  for  -2.5°  declination  are  obtained.  Conversely,  where  the 
declination  of  the  field  line  is  eastward  in  the  southern  hemisphere,  a 
westward-blowing,  early  morning  neutral  wind  transports  ionization  downward 
into  a  higher  loss  rate  region,  decreasing  the  peak  electron  density  such  as 
slnwn  for  18°  declination.  During  nighttime  hours  when  the  zonal  wind  blows 
eastward,  the  effer’-s  of  field-line  declination  on  peak  electron  density  are 
reversed.  In  the  southern  hemisphere,  an  eastward-blowing  wind  lowers  the  F 
region  if  the  declination  is  westward,  thus  decreasing  foF’  while  same 
eastward  wind  raises  the  F  layer  where  declination  is  eastward,  causing  a 


~o  DECL  =  -35.4°,  DIP  =  -69.4°,  Bt.6°  E.  6E0G  LONG,  40.6  S.  GEOG  LAT 


jguro  *>.  Diurnal  behavior  of  foF2  during  September  1978  deduced  using  the  continuity  equation 
with  the  same  neutral  wind  model  at  three  locations  in  the  southern  hemisphere 
having  essentially  the  same  dip  but  different  declinations. 


slower  decay  in  foF2  with  time.  This  can  be  seen  by  comparing  the  slopes  of 
the  curves  shown  in  Figure  5  between  the  hours  1800  and  2200  local  time. 

The  foregoing  illustrates  how  the  theoretical  model  provides  a  convenient 
and  valid  method  to  obtain  physically  realistic  values  of  foF2  at  locations 
that  are  generally  inaccessible  to  routine  observation. 

4.  MATHEMATICAL  BASIS  FOR  THE 

GENERATION  OF  THE  foF2  COEFFICIENTS 

The  mathematical  basis  for  generating  the  numerical  coefficients  that 
represent  the  global  variation  of  foF2  has  been  provided  by  Jones  and  Gal  let 
(1962).  The  application  of  the  Jones  and  Gal  let  (1962)  method  to  ionospheric 
data  has  been  discussed  by  numerous  authors,  most  recently  by  Rush  et  al . 
(1982).  In  this  section,  a  brief  description  of  the  methodology  used  in 
developing  the  numerical  coefficients  is  given.  More  detailed  discussion  can 
be  obtained  from  Jones  et  al.  (1969)  and  Jones  and  Obitts  (1970). 

In  choosing  a  functional  representation  for  foF2,  consideration  must  be 
given  to  the  nature  of  the  data  itself  and  the  large-scale  variations  that 
must  be  represented.  Because  the  F2-region  critical  frequency  displays  a 
three  dimensional  (time,  geographic  latitude,  geographic  longitude)  variation, 
a  simple  set  of  functions  will  not  represent  the  behavior  of  foF2.  Also,  the 
geographic  distribution  of  the  data  used  to  derive  the  numerical  coefficients 
presents  problems  in  preserving  the  stability  of  the  final  result  in  areas 
where  few  data  points  (or  observations)  are  available.  Ideally  equally  spaced 
data  distributed  uniformly  over  the  globe  and  in  time  are  desired  in  order  to 
obtain  the  final  coefficients. 

The  usual  method  of  representing  global  variations  of  geophysical  data  is 
to  use  spherical  harmonic  analysis.  This  approach  provides  functions  that  are 
periodic  in  longitude  and  constant  at  the  poles.  The  classical  surface 
spherical  harmonics  that  satisfy  Laplace's  equations,  however,  are  not  ortho¬ 
gonal  with  respect  to  the  locations  for  which  ionosphere  data  are  available. 
Therefore,  two  sets  of  functions  were  chosen,  A  set  of  orthogonal  polynomials 
was  chosen  to  represent  geographic  variations.  The  diurnal  variations  then 
need  only  be  represented  by  a  periodic  function,  and  a  Fourier  series  was 
chosen  to  provide  this  representation.  Advantages  of  using  Fourier  series  are 
that  adjusting  the  data  for  differences  in  recording  times  is  vgreatly 
simplified,  the  noise  inherent  in  the  observations  is  more  easily  separated 


from  the  true  diurnal  variation  of  foF2,  and  the  resultant  accuracy  of  the 
Fourier  fit  to  the  data  can  be  readily  assessed  using  least-squares 
techniques. 

The  application  of  Fourier  analysis  to  represent  the  diurnal  variation  of 
foF2  using  k  harmonics  yields  an  expression  of  the  form 


where 


fk(x)  =  aQ  +  £  [2j  cos  jx  +  b .  sin  j 

j — ^ 


jx] 


(2) 


-(5s)  -  •*. 


i  =  0,  1,  .  .  .  N, 


1  N 

*o  =TT  ^  yr 


(3) 


3j  =  Tf  £  y1  C0S  jxi 


jnd  bj  =4  stn  Jxt 

for  all  values  of  k  where  (2k  +  1  £  N).  The  value  of  Yk(x)  is  a  best 
approximation  in  the  sense  of  least  squares.  In  the  current  application,  the 
24  hourly  measurements  of  foF2  are  available  usually  in  local  standard  time 
and  therefore  must  be  translated  to  local  time.  Eight  harmonics  were  used  in 
the  Fourier  analysis  which  agrees  with  previous  studies  (Rush  et  al .  (1982). 
This  tends  to  smooth  out  noise  and  at  the  same  time  gives  a  good  representa¬ 
tion  of  the  diurnal  variation  of  foF2.  After  fitting  the  diurnal  data  with 
the  Fourier  series,  the  Fourier  coefficients  are  shifted  in  time  so  that  the 
data  at  each  station  can  be  represented  in  local  time.  The  Fourier  coeffi¬ 
cients  representing  the  diurnal  variation  of  the  foF2  data  are  then  fit  with 
respect  to  the  geographic  location  at  which  the  data  were  taken. 

The  set  of  functions  chosen  to  represent  the  geographic  variation  of  the 
Fourier  coefficients  is  given  in  Table  4.  This  set  of  geographic  functions 


Table  4.  Set  of  Functions  Chosen  to  Represent 
the  Geographic  Variation  of  the 
Fourier  Coefficients 


Pi 

?l  +  1 

Pi  +  2 
Pi  +  3 
px  +  4 


cosA  cose 
cosA  sine 

sinX  cosA  cose 

sinX  cosA  sine 


Pz-  i 

P2 

p2  +  1 
p2  +  2 

P2  +  3 
P2+  4 


p3  -  1 

P3 


sinaX  cosA  cose 
sinaX  cosA  sine 
Cos2A  cos2e 
cos2A  sin2e 
sinX  cos2A  cos20 

sinX  cos2A  sin20  b  =  P,3  Pi  -  1 

2 

« 

sinbX  cos2X  cos2e 
sinbX  cos2A  sin26 


where  A  =  Geographic  latitude, 

8  =  Geographic  longitude, 
sinX  =  I/(l2  *  cosA)1^2, 
and  I  =  magnetic  inclination  (1975  epoch). 
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G^ ( X,  0)  given  in  Table  4  is  analogous  in  some  respects  to  spherical  harmonics. 
The  main  latitudinal  variation  terms  correspond  to  zonal  harmonics,  and  the 
mixed  latitudinal  and  longitudinal  terms  Involving  cos  &  or  sin  9  (along  with 
the  terms  involving  cos  20  and  sin  20)  correspond  to  the  tesseral  harmonics  of 
the  first  and  second  order,  respectively.  The  series  of  functions  G^(A,  0) 
are  periodic  in  longitude,  constant  at  the  poles,  and  are  weighted  according 
to  latitude  by  the  functions  cosX,  cos^X',  .... 

In  terms  of  orthogonal  polynomials,  the  desired  geographic  functions  may 
be  written  simply  as  (Rush  et  al.,  1982) 

\(x)  =  do  PoW  =  dl  P1  M  *  +  dk  Pk  <x>  (4) 

where  the  coefficients  in  the  orthonogal  system  are 
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The  numbers  Cj  are  given  by 


I  [Pi  (x,)]2  ,  and 
>1  J  1 


the  coefficients  in  the  orthonormal  system  are 

V4j*t  <7> 

Each  term  djPj(x)  in  the  series  (8)  is  independent  of  the  other  terms.  As  a 
consequence,  the  degree  of  the  approximating  polynomial  may  be  increased  with¬ 
out  affecting  the  terms  already  obtained  and,  in  some  cases,  tests  for  statis¬ 
tical  significance  may  be  applied  to  each  term. 

The  accuracy  of  the  coefficients  determined  from  both  the  diurnal  and 
geographic  analysis  was  computed  by  use  of  a  least-squares  method  described  by 
Rush  et  al .  (1983).  The  root-mean-square  (rms)  error  for  the  fit  to  the  foF2 
data  (observation,  contour-derived  low-latitude  data,  and  theoretically- 
derived  higher  latitude  data)  for  each  month  for  the  solar  minimum  (1975/1976) 
and  solar  maximum  (1978/1979)  time  periods  is  given  in  Table  5.  Also  shown  in 
the  table  is  the  rms  error  by  Jones  et  al.  (1969)  in  deriving  the  CCIR  numeri¬ 
cal  coefficients  fo*'  the  F2-region  critical  frequency.  The  rms  error  for  the 
new  coefficients  is  on  the  order  of  0.4  to  0.5  MHz  for  the  solar  minimum 
representation  and  between  0.5  and  0.7  MH2  for  the  solar  maximum  representa¬ 
tion,  It  can  be  seen  that  for  the  new  representation,  the  rms  errors  during 


Table  5.  RMS 
to 

Errors  (MHz)  for  the  Fit  to  the  Data 
Generate  foF2  Coefficients 

Used 

Month 

New  Coefficient 

Solar  Minimum  Solar  Maximum 

Jones  et  al 
Solar  Minimum 

.  (1969) 

Solar  Maximum 

Jan 

0.513 

0.558 

.406 

.662  j 

Feb 

0.375 

0.604 

.360 

.607  ! 

Mar 

0.428 

0.631 

.378 

.619  ' 

Apr 

0.425 

0.633 

.379 

.600  ! 

May 

0.310 

0.622 

.344 

.585 

Jun 

0.334 

0.481 

.352 

.583 

Jul 

0.337 

0.474 

.332 

.581 

Aug 

0.392 

0.516 

.353 

.553 

Sep 

0.449 

0.619 

.404 

.589 

Oct 

0.452 

0.632 

.402 

.650 

Nov 

0.400 

0.703 

.415 

.699 

Dec 

0.379 

0.560 

.418 

.696 

the  solar  minimum  period  are  generally  one-half  to  two-thirds  of  those  at 
solar  maximum.  Also,  the  errors  tend  to  be  smaller  during  the  northern  hemi¬ 
sphere  summer  months  (May  through  August)  than  for  the  other  months  of  the 
year. 

The  rms  errors  are  comparable  to  those  obtained  from  the  Jones  et  al . 
(1969)  analysis  and  are  probably  close  to  the  noise  level  of  the  data  used  in 
the  generation  of  the  coefficients.  The  errors  could  be  reduced  somewhat  by 
including  higher  order  harmonics  into  the  spherical  harmonic  analysis.  The 
use  of  higher  order  harmonics  to  represent  the  observed  foF2  data  could, 
however,  give  rise  to  instabilities  in  the  representation  of  foF2  in  data- 
sparse  regions.  Also,  the  coefficients  generated  would  not  be  of  the  same 
format  as  the  CCIR  coefficients  (CCIR,  1982b).  One  of  the  desired 
accomplishments  at  the  onset  of  this  study  was  to  produce  a  more  realistic  set 
of  numerical  coefficients  to  represent  foF2  on  a  global  scale  that  could  be 
directly  substituted  for  the  CCIR  coefficients.  Thus,  the  decision  to 
maintain  the  same  number  of  harmonics  as  given  by  the  CCIR  was  take*-. 

The  entire  set  of  coefficients  that  represent  the  global  variation  of 
foF2  is  listed  in  Appendix  il  for  each  month  of  July  1975  to  June  1976  and  July 
1978  to  June  1979. 


24 


5.  RESULTS  AND  DISCUSSION 

The  numerical  coefficients  that  represent  the  global  variation  of  foF2 
were  derived  from  the  data  described  in  Sections  2  and  3  using  the  procedures 
outlined  in  Section  4.  For  each  month  between  July  1975  and  June  1976  and 
between  July  1978  and  June  1979,  a  set  of  coefficients  has  been  derived  that 
is  based  on  observations  of  foF2  available  from  the  worldwide  network  of  ver¬ 
tical  incidence  ionosonde  stations.  In  this  study,  these  observations  have 
been  supplemented  by  data  obtained  from  a  detailed  constant-local -time- 
analysis  of  foF2  observations  at  low  magnetic  latitude  and  by  data  obtained 
from  the  time-dependent  continuity  equation  at  the  higher  latitudes.  Typi¬ 
cally,  for  each  month  about  20  locations  with  actual  observations  and  95  loca¬ 
tions  of  contour-derived  foF2  values  were  used  for  the  low-latitude  data 
base.  For  the  higher  latitudes  for  each  month,  about  75  locations  of  actual 
observations  and  90  locations  of  theoretically  derived  data  formed  the 
re  uired  data  base.  If  a  value  of  foF2  determined  from  the  continuity  equa¬ 
tion  was  below  0.7  MHz,  the  value  was  set  equal  to  0.7  MHz.  This  was  done  in 
order  to  assure  that  the  values  of  foF2  subsequently  determined  from  the 
numerical  coefficients  did  not  yield  F2-region  critical  frequencies  that  were 
near  zero  or  negative. 

An  example  of  the  global  variation  of  foF2  obtained  from  the  coefficients 
derived  in  this  study  is  shown  in  Figure  6.  Monthly  median  contours  of  foF2 
(in  MHz)  are  illustrated  for  January  1976  at  1200  hrs  universal  time  (>JT). 

The  figure  provides  a  striking  example  of  the  high  latitude  trough  of  ioni¬ 
zation  typical  of  the  winter  nighttime  hours  (note  the  contours  between  90° 
and  270°  longitude,  and  poleward  of  40°  north  latitude).  Also  shown  for  the 
late  morning  and  afternoon  hours  (between  315°  and  120’’  longitude)  at  latitudes 
between  30°  north  and  south  is  the  equatorial  anomaly.  This  region  is  charac¬ 
terized  by  ionization  maxima  bounding  a  relative  minima  of  ionization  above 
the  magnetic  equator.  The  summertime  hemisphere,  particularly  for  middle 
southern  hemispheric  latitudes,  displays  values  of  foF2  that  vary  little  with 
longitude  (and  local  time).  This,  too,  is  typical  of  the  F2  region  morphol¬ 
ogy.  The  major  ionization  features  depicted  on  the  figure  are  aligned  with 
magnetic  latitude  which  explains  why,  in  the  geographic  representation  of  the 
figure,  the  ionization  features  tend  to  be  slanted  with  respect  to  geographic 
latitude. 
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Contour  map  of  the  global  representation  of  the  median  value  of  foF2  for  January 
1976  at  1200  hrc  UT  derived  from  the  new  coefficients. 


Figures  7,  8,  and  9  show  contours  of  foF2  for  1200  hrs  for  April  1976, 
July  1975,  and  October  1975,  The  months  of  January,  April,  July,  and  October 
were  chosen  for  illustration  since  these  months  are  the  ones  indicated  by  the 
CCIR  as  being  the  most  representative  of  the  seasonal  variations  of  the  iono¬ 
sphere  (CCIR,  1983).  These  figures,  taken  together  with  Figure  6,  provide  an 
illustration  of  the  seasonal  variations  of  foF2  during  a  period  of  minimum 
solar  activity.  It  is  not  the  purpose  of  this  report  to  describe  the  morphol¬ 
ogy  of  the  F2  region  that  is  given  by  the  numerical  representation  of  foF2. 
Suffice  it  to  say  that  the  general  features  of  the  global  variation  of  foF2 
given  in  Figures  6  through  10  are  in  total  accord  with  the  known  morphology  of 
foF2  that  has  emerged  as  a  result  of  over  50  years  of  ionospheric  research. 

In  Appendix  C,  global  maps  of  foF2  are  given  for  0600  and  1800  hrs  UT  for  the 
four  seasonal  months  for  both  solar  minimum  (1975/1976)  and  solar  maximum 
(1978/1979).  The  maps  given  in  Appendix  C  illustrate  the  seasonal  and  solar 
cycle  dependencies  observed  in  foF2. 

Figures  10  and  11  illustrate  the  global  variation  of  foF2  for  January  and 
April  1976  at  l?on  hrs  UT  derived  from  the  numerical  coefficients  given  by 
Jones  et  al.  (1969).  These  coefficients,  referred  to  as  the  CCIR  coeffi¬ 
cients,  provide  a  representation  of  foF2  that  is  based  on  data  observed 
between  1954  and  1958.  More  Importantly,  however,  is  the  fact,  that  for 
regions  of  the  earth  where  data  were  unavailable,  simple  interpolation  pro¬ 
cedures  were  used  to  obtain  foF2  values.  Comparing  Figures  10  and  11  with  6 
and  7  shows  that  the  high-latitude  trough  in  Figures  6  and  7  is  much  better 
defined  than  in  Figures  10  and  11.  Also,  comparing  Figures  6  and  7  with 
Figures  10  and  11,  it  is  seen  that  the  form  of  the  foF2  contours  in  the 
southern  hemisphere  differs  substantially.  This  is  due  to  the  fact  that  the 
data  used  to  generate  Figures  6  and  7  were  derived  from  the  continuity  equa¬ 
tion  and  reflect  the  dependence  of  foF2  on  both  magnetic  declination  and 
inclination.  The  data  used  to  generate  Figures  10  and  11,  on  the  other  hand, 
only  reflect  the  magnetic  inclination  (or  dip)  control  of  the  F2  region. 

While  comparisons  such  as  those  discussed  above  are  useful  for  illustra¬ 
tive  purposes,  more  detailed  comparisons  between  actual  observations  and  foF2 
derived  from  the  new  coefficients  and  the  CCIR  coefficients  are  necessary  to 
determine  what,  if  ary,  improvement  in  the  global  representation  of  foF2  is 
afforded  by  the  new  coefficients.  Befote  discussing  improvements,  however,  it 
is  worthwhile  to  consider  the  differences  in  the  values  of  foF2  determined 
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Contour  map  of  the  global  representation  of  the  median  value  of  foF2  for  April  1976 
at  1200  hrs  UT  derived  from  the  new  coefficients. 
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Contour  map  of  the  global  representation  of  the  median  value  of  foF2  for  July  1975 
at  1200  hrs  UT  derived  from  the  new  coefficients. 


Contour  map  of  the  global  representation  of  the  median  value  of  foF2  for  October 
1975  at  1200  hrs  UT  derived  from  the  new  coefficients. 


from  the  two  sets  of  coefficients  themselves.  Figure  12  shows  the  diurnal 
variation  of  foF2  for  March  1979  that  has  been  derived  from  the  new 
coefficients  and  that  derived  from  the  existing  CCIR  coefficients.  Values  of 
foF2  are  shown  for  three  distinct  ranges  of  declination  (eastward,  westward, 
and  approximately  zero  declination)  and  for  five  different  latitude  intervals 
(northern  middle,  northern  low,  equatorial,  southern  low,  and  southern  middle 
latitude).  The  latitudinal  intervals  are  designated  according  to  magnetic 
latitude.  The  coordinates  of  the  15  locations  given  in  the  figure  are  listed 
in  Table  6.  For  those  locations  where  observations  of  foF2  were  available  in 
the  generation  of  the  numerical  coefficients,  the  station  indicator  begins 
with  a  numeral  and  the  station  name  is  listed  in  Table  6. 

The  differences  between  foF2  derived  from  the  two  sets  of  coefficients 
are  rather  striking  at  many  of  the  locations.  This  is  particularly  the  case 
for  those  locations  where  observations  of  foF2  were  not  available  (i.e.,  the 
locations  beginning  with  a  letter).  It  is  at  these  locations  that  the  biggest 
differences  in  the  values  of  foF2  derived  from  the  two  sets  of  coefficients 
would  be  expected.  The  reason  for  this  is  that,  at  unobserved  locations,  the 
values  of  foF2  derived  from  the  numerical  coefficients  is  dependent  upon  the 
procedure  used  to  fill  in  data  at  the  locations  in  question.  As  was  mentioned 
earlier,  the  CCIR  coefficients  were  generated  from  data  that  were  simply 
interpolated  (or  extrapolated)  taking  account  of  only  the  magnetic  inclination 
control  of  the  F2  region.  The  new  coefficients,  on  the  otner  hand,  were  gen¬ 
erated  from  data  that  accounted  for  the  magnetic  inclination  and  magnetic 
declination  control  of  the  F2  region. 

The  results  shown  in  Figure  12  illustrate  quite  readily  the  impact  of 
neutral  winds  on  the  diurnal  variation  of  foF2  determined  from  the  new  coeffi¬ 
cients  at  all  locations.  This  is  particularly  noticeable  for  the  northern  and 
southern  middle-latitude  locations.  For  the  051,  945,  and  336  location  group¬ 
ing,  the  critical  frequency  during  the  day  hours  is  highest  for  the  eastward 
declination  and  during  the  night  foF2  is  highest  for  the  westward  declination. 
The  critical  frequency  for  the  85!,  B6M,  and  Y5!  location  grouping  shows 
exactly  the  opposite  behavior:  highest  values  of  foF2  during  the  day  are 
found  at  the  westward  doc ii nation  and  during  the  night  at  the  eastward  declin¬ 
ation.  This  behavior  results  from  the  effect  of  the  zonal  component  of  the 
neutral  air-wind  on  the  variation  of  foF2.  It  can  be  seen  from  Figure  A. 3 
that  the  zonal  wind  is  directed  eastward  in  the  evening  and  westward  in  the 
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Table  6.  Coordinates  of  Locations  Used  to  Illustrate 
Differences  in  foF2  Derived  from  the  CCIR 
Coefficients  and  the  New  Coefficients 


Location 

Indicator 

Geog. 

Latitude 

(°NL 

Geog. 

Longitude 

<°E> 

Magnetic 

Inclination 

HO 

Magnetic 

Declination 

(°E) 

Location 

Name 

051 

51.5 

359.4 

66.8 

-  7.0 

SI ough 

•  v; 

945 

45.0 

285.0 

75.0 

0.0 

Ottawa 

ri  •  .  , 

836 

35.6 

239.4 

60.7 

15.0 

Point  Arguello 

‘ 

L15 

5.0 

310.0 

25.1 

-16.0 

L07 

0.0 

290.0 

23.9 

-  2.0 

L06 

0.0 

240.0 

9.9 

8.0 

L2N 

-10.0 

310.0 

-6.0 

-15.0 

91K 

-12.0 

285.0 

1.0 

2.0 

Huancayo 

L05 

0.0 

170.0 

3.5 

10.0 

L3L 

-15.0 

340.0 

-47.2 

-22.0 

,  -/* .  •. 

L6K 

-30.0 

300.0 

-30.3 

-  2.0 

L5J 

-25.0 

250.0 

-33.6 

15.0 

.  * . 

** » *  •  • . 

B5! 

-49.6 

344.4 

-57.8 

-22.0 

w 

^  .  1,  ».  . 

B6M 

-64.3 

316.5 

-58.3 

5.0 

• 

• 

Y5! 

-50.2 

236.2 

-62.1 

27.0 

fo  F2  (MHz) 


Figure 


East  declination  Zero  declination  West  declination 


March  1979 


.  Diurnal  variation  of  the  median  value  of  foF2  during  March  1979 

determined  from  the  new  coefficients  and  the  CC I R  coefficients  for  the 
northern  middle,  northern  low,  equatorial,  southern  low,  and  southern 
middle  latitudes  at  locations  with  westward,  eastward,  and  zero  declin 
ations. 
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early  morning  hours  in  both  hemispheres.  The  zonal  wind  is  directed  westward 
during  the  day  in  the  southern  hemisphere  and  is  zero  during  the  day  in  the 
northern  hemisphere.  In  the  northern  hemisphere,  the  westward  wind  during  the 
early  morning  hours  drives  ionization  up  the  magnetic  field  lines  where  the 
declination  is  eastward  and  down  the  field  lines  where  the  declination  is 
westward  (see  Figure  3).  The  ionization  that  is  moved  upward  is  move':’  into 
regions  of  lower  loss  rates.  The  values  of  foF2  are  therefore  larger  at  these 
locations  than  at  locations  where  the  ionization  is  moved  downward  into  higher 
loss  regions.  In  the  southern  hemisphere,  recalling  Figure  3,  the  opposite 
effect  is  noted.  Similar  arguments  can  be  made  concerning  the  effect  of  the 
zonal  wind  on  foF2  at  the  other  local  times. 

Figure  12  also  provides  evidence  that  neutral-air  winds  affect  the 
diurnal  behavior  of  the  ionization  distribution  at  the  low  latitudes  and  even 
near  the  equator.  For  example,  fof2  at  the  northern  low-latitude  locations 
during  the  early  morning  hours  is  noticeably  lower  at  the  westward  declination 
than  at  the  zero  declination  location.  This  agrees  with  the  expected  impact 
of  the  zonal  wind  in  the  northern  hemisphere  which  is  directed  westward  during 
the  early  morning  hours.  Such  effects  observed  at  the  low  latitudes  are  par¬ 
ticularly  encouraging  because  the  data  used  to  supplement  the  low-latitude 
observations  were  not  generated  using  the  time-dependent  continuity  equation. 
The  continuity  equation  yields  values  of  foF2  that  must  agree  with  the  effect 
of  neutral-air  winds  on  she  ionization  distribution.  The  low-latitude  coeffi¬ 
cients  were  generated  from  data  that  were  either  observed  or  determined  from 
the  low-latitude  local  time  analysis.  That  the  resultant  values  of  foF2  agree 
with  theoretical  expectation  provides  further  evidence  of  the  validity  of  the 
data  obtained  from  the  local  time  analysis. 

Figure  13  provides  results  for  the  solar  minimum  period  of  March  1976 
similar  to  those  shown  in  Figure  12.  The  results  displayed  in  Figure  13  are 
also  in  agreement  with  theoretical  expectation.  Figure  A. 3  and  Figure  3  can 
be  used  to  verify  that  this  is  indeed  the  case. 

In  order  to  determine  what  improvement  in  foF2  is  gained  by  using  numer¬ 
ical  coefficients  that  are  consistent  with  theoretical  expectation,  differ¬ 
ences  between  observations  of  foF2  and  foF2  deduced  from  both  sets  of  coeffi¬ 
cients  (the  new  coefficients  and  the  CCIR  coefficients)  must  be  compared.  It 
is  extremely  difficult  to  obtain  data  that  can  be  used  to  objectively  compare 
the  values  of  fof2  with  observations,  particularly  in  regions  of  the  globe 
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new  coefficients  and  the  CCIR  coefficients  for  the  northern 
equatorial ,  southern  low,  and  southern  middle  latitudes  at 
eastward,  and  zero  declinations. 


where  it  is  expected  that  the  new  coefficients  will  show  the  greatest 
improvement,  such  as  over  the  oceans.  The  reasons  for  this  are  rather 
obvious.  Any  location  that  is  far  removed  from  other  locations  and  for  which 
data  were  available  was  used  in  the  generation  of  the  coefficients.  To  do 
otherwise  would  result  in  coefficients,  and  hence  values  of  foF2,  that  could 
yield  needlessly  large  errors  at  the  very  locations  where  improvement  in  the 
specification  of  foF2  is  needed  most.  One  possible  approach  to  overcome  this 
dilemma  is  to  use  the  numerical  coefficients  to  predict  values  of  foF2  and 
then  to  compare  the  predicted  values  with  observations. 

The  accuracy  of  the  predicted  values  of  foF2  depends  not  only  on  the 
values  of  foF2  specified  for  each  month  by  the  coefficients  but  also  on  the 
relationship  between  solar  activity  and  foF2.  The  same  relationship  between 
foF2  and  solar  activity  was  assumed  for  values  determined  from  both  sets  of 
coefficients.  This  relationship  is  a  linear  one  between  foF2  and  the  12-month 
average  sunspot  number,  Rl2>  given  by 


where  Ry  =  value  of  12-month  smoothed  sunspot  number  for  solar  maximum 
conditions , 

R  =  value  of  12-month  smoothed  sunspot  number  for  solar  minimum 
conditions, 

=  value  of  foF2  at  R^, 

Fu  =  value  of  foF2  at  Ry, 

and  R^2  =  value  of  12-month  smoothed  sunspot  number  for  month 
predictions  are  being  made. 

It  was  further  assumed  that  any  errors  in  foF2  resulting  from  an  error  in  this 
linear  relationship  will  effect  both  sets  of  coefficients  equally.  The  valid¬ 
ity  of  this  assumption  is  questionable,  and  further  study  of  this  is  certainly 
warranted. 

Figures  14,  IS,  and  16  are  maps  of  foF2  determined  from  the  -new  coeffi¬ 
cients,  showing  contours  of  the  critical  frequency  at  0000  hrs  UT  in  July  for 
solar  conditions  appropriate  to  sunspot  numbers  of  15  (Figure  14),  55 
(Figure  15),  and  95  (Figure  16).  The  results  for  Figure  14  are  valid  for  July 
1975  and  those  of  Figure  16  are  valid  for  July  1973.  The  values  of  fo?2  given 


33 


Contour  map  of  the  global  representation  of  the  med 
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in  Figure  15  were  determined  by  linearly  interpolating  between  the  coef¬ 
ficients  for  July  1975  and  July  1978.  It  is  seen  from  the  three  figures  that 
the  general  pattern  of  the  variations  in  foF2  remains  the  same  as  the  level  of 
solar  activity  increases.  The  values  of  the  contours  tend  to  increase  as 
solar  activity  increases.  Values  of  foF2  determined  from  representations  such 
as  illustrated  in  Figure  15  provided  the  basis  for  the  comparisons  between 
foF2  given  by  both  sets  of  coefficients  and  by  observations  of  foF2. 

Ideally,  in  order  to  compare  the  accuracy  of  the  predictions  of  foF2 
using  the  new  coefficients  and  using  the  current  CCIR  coefficients,  predic¬ 
tions  should  be  made  for  time  periods  and  for  locations  that  were  not  used  in 
the  development  of  either  set  of  coefficients.  The  new  coefficients  presented 
here  were  developed  using  data  observed  during  and  appropriate  to,  the  solar 
minimum  year  of  July  1975  through  June  1976  and  the  solar  maximum  year  of  July 
1978  through  June  1979.  The  CCIR  coefficients  have  been  developed  using  data 
valid  for  the  solar  minimum  period  of  1954  and  the  solar  maximum  period  of 
1958.  In  order  to  be  as  objective  as  possible,  comparisons  of  the  prediction 
accuracy  using  both  sets  of  coefficients  were  made  for  years  spanning  the 
early  1960's  through  1972. 

Table  7  shows  the  rms  prediction  errors  in  foF2  determined  at  Freiburg 
(48. 1°N,  7.6°E)  for  the  months  of  January,  April,  July,  and  October  for  the 
years  1964  through  1971.  Data  from  Freiburg  were  not  used  in  the  generation 
of  either  set  of  coefficients.  It  can  be  seen  that  for  the  months  studied,  in 
20  out  of  32  cases,  the  use  of  the  new  coefficients  to  predict  foF2  at 
Freiburg  resulted  in  smaller  prediction  errors  than  using  the  CCIR  coeffi¬ 
cients.  Also  worth  noting  is  that  the  new  coefficients  tend  to  have  lower 
prediction  errors  than  the  CCIR  coefficients  during  t  ie  solar  minimum  years  of 
1964,  1965,  1970,  and  1971  (13  out  of  the  16  months)  bit  not  during  tHe  solar 
maximum  years  of  1965,  1967,  1968,  and  1969  (7  out  of  ti.  16  mointh*-'^ 

Comparisons  for  other  locations  show  the  1  use  of  the  new  -^efficients  do 
not  always  lead  to  improved  predictions  of  foF2.  Table  8  gives  rms  prediction 
errors  of  foF2  using  both  the  new  coefficients  and  the  CCIR  coefficients  for 
Grand  Bahama,  Singapore,  Cocos  Island,  and  Hon  ^ong.  Comparisons  were  made 
for  the  same  four  months  as  given  in  Table  7  t  for  the  solar  minimum  year  of 
1964  and  the  solar  maximum  year  of  1969.  The  prediction  errors  indicate  that 
toe  new  coefficients  do  not  improve  the  prediction  of  to?2,  particularly 
dvring  solar  maximum  conditions  (1969).  It  sh mid  be  rioted, 


Table  7.  RMS  Prediction  Errors  (MHZ)  in  foF2  for  Freiburg 


Coefficient  Set 

JAN 

APR 

JUL 

OCT 

1964 

New 

.54 

.34 

.21 

.32 

CCIR 

.35 

.45 

.23 

.45 

1965 

New 

.65 

.47 

.23 

.27 

CCIR 

.29 

.64 

.28 

.33 

1966 

New 

.56 

.39 

.26 

.92 

CCIR 

.40 

.49 

.63 

.41 

1967 

New 

.42 

1.05 

.39 

.49 

CCIR 

.33 

1.23 

.48 

.65 

1968 

New 

.49 

.32 

.67 

.78 

CCIR 

.94 

.31 

.29 

.20 

1969 

New 

.79 

.52 

.50 

.95 

CCIR 

.60 

.69 

.45 

.32 

1970 

New 

.52 

.99 

.24 

.45 

CCIR 

.77 

1.18 

.68 

.61 

1971 

New 

.48 

.31 

.27 

1.09 

CCIR 

.61 

.41 

.42 

.55 

•n 
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Table  8.  RMS  Prediction  Errors  (MHZ)  in  foF2  at  Selected 

Locations  for  1964  and  1969  Solar  Activity  Conditions 


Coefficient 


Set 

JAN 

APR 

JUL 

OCT 

1964  - 

New 

.90 

.43 

.39 

,53 

CCIR 

.48 

.32 

.30 

.37 

GRAND  BAHAMA 

1969 

New 

1.40 

.48 

.66 

1.05 

CCIR 

.99 

.50 

.41 

.65 

1964 

New 

.53 

.50 

OJ 

• 

.51 

CCIR 

.54 

.37 

.45 

.69 

SINGAPORE 

1969 

New 

1.31 

.62 

1.12 

1.19 

however,  that  data  for  three  of  the  four  locations  given  in  Table  8--Grand 
Bahama,  Singapore,  and  Cocos  Island--were  used  in  the  development  of  the  CCIR 
coefficients  but  not  in  the  development  of  the  new  coefficients.  This  may 
explain,  in  part  at  least,  the  better  agreement  using  the  CCIR  coefficients. 

Table  9  shows  the  results  of  a  comparison  of  predicted  foF2  values  for 
each  month  of  the  year  of  1970  at  seven  locations  on  the  earth  whose  data  were 
used  in  the  generation  of  both  sets  of  coefficients.  The  locations  shown  in 
the  table  are  Slough,  England;  Tokyo,  Japan;  Maui,  Hawaii;  Huancayo,  Peru; 
Tahiti,  Pacific  Ocean;  Canberra,  Australia;  and  Christchurch,  New  Zealand.  It 
can  be  seen  that  at  Slough,  Huancayo,  Tahiti,  and  Canberra,  the  new  coeffi¬ 
cients  yield  values  of  foF2  that  are  generally  closer  to  the  observed  values 
than  do  the  CCIR  coefficients.  On  the  other  hand,  at  Tokyo  and  Maui,  the  CCIR 
coefficients  yield  predictions  that  are  closer  to  the  observations  than  are 
the  new  coefficients  for  most  months  of  the  year. 

It  is  highly  likely  that  the  results  presented  above  are  colored  by 
errors  in  the  solar  activity  dependence  assumed  for  the  new  coefficients  in 
equation  (8).  The  apparent  lack  of  overall  improvement  in  the  prediction  of 
foF2  using  the  new  coefficients  may  be  due  more  to  deficiency  in  the  solar 
activity  interpolation  procedure  than  to  a  deficiency  in  the  coefficients 
themselves.  This  is  based  on  the  fact  that  for  solar  minimum  conditions  the 
new  coefficients  tend  to  yield  predicted  values  of  foF2  that  are  either  as 
good  as  or  better  than  those  given  by  the  CCIR  coefficients.  The  solar  activ¬ 
ity  for  the  months  used  to  determine  the  new  coefficient  for  solar  minimum 
was,  in  fact,  indicative  of  minimum  solar  activity  (see  Table  1).  Thus, 
predictions  of  foF2  at  solar  minimum  are  based  on  coefficients  that  are 
realistic  representations  of  solar  minimum  conditions.  The  same  cannot  be 
said,  however,  for  the  new  coefficients  for  solar  maximum.  The  values  of  the 
sunspot  numbers  used  to  determine  the  new  coefficients  varied  from  95  in  July 
1978  to  over  150  in  June  1979.  While  values  of  sunspot  numbers  of  140  to  150 
are  indicative  of  solar  maximum  conditions,  values  of  95  to  120  may  not  be. 

The  values  of  foF2  determined  from  the  new  coefficients  for  solar  maximum 
conditions  may  be  too  low  because  they  are  based  on  data  that  may  not  be 
totally  representative  of  solar  maximum  conditions. 

The  results  presented  in  Table  7  clearly  support  the  above  argument.  It 
would  be  expected  that  the  predictions  of  foF2  for  the  month  of  April  that  are 
based  on  the  new  coefficients  would  be  more  accurate  than  the  CCIR 
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Table  9.  RMS  Prediction  Errors  (MHZ)  in  foF2  for  the 
Locations  Given  for  Each  Month  of  1970 
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coefficients  because  the  data  used  to  derive  the  new  soiar  maximum  coeffi¬ 
cients  for  April  were  applicable  to  a  sunspot  number  (14t)  that  is  truly 
representative  of  solar  maximum  conditions.  Of  the  results  presented  for 
April  in  Table  7,  seven  out  of  eight  cases  show  that  for'2  determined  from  the 
new  coefficients  is  better  than  foF2  determined  from  the  CCIR  coefficients. 
Even  in  Table  8,  the  results  given  for  the  month  of  April  show  that  the  CCIR 
coefficients  do  not  yield  values  of  foF2  that  are  vastly  better  than  those 
determined  from  the  new  coefficients.  This  is  particularly  noteworthy 
because,  as  mentioned  previously,  observations  at  Grand  Bahama,  Singapore,  and 
Cocos  Island  were  used  in  the  generation  of  the  CCIR  coefficients  but  not  in 
the  generation  of  the  new  coefficients. 

The  results  presented  above  indicate  that  the  values  of  foF2  obtained 
from  the  new  coefficients  appear  to  be  better  than  those  obtained  from  the 
CCIR  coefficients  if  proper  account  can  be  taken  of  the  solar  cycle  dependence 
of  the  coefficients.  This  result  holds  true  for  regions  of  the  globe  where 
observations  of  foF2  were  generally  available  to  determine  the  numerical 
coefficients.  At  locations  where  data  are  not  available,  it  is  difficult  to 
quantify  the  improvement  afforded  by  the  new  coefficients.  However,  Rush  et 
al.  (1983)  have  shown  that,  for  the  southern  hemisphere,  values  of  foF2 
determined  from  the  new  coefficients  are  generally  closer  to  those  obtained 
from  the  Japanese  satellite  ISS-b  than  are  the  values  of  foF2  deduced  from  the 
CCIR  coefficients  during  December  1978  (a  period  of  solar  maximum). 

Figure  17  (which  is  Figure  11  in  Rush  et  al . ,  1982)  shows  three  longitu¬ 
dinal  distributions  of  t‘oF2  at  three  latitudes  (3C°S,  40°5,  and  50°S)  for 
December  1978,  1000  hrs  UT.  Shown  is  foF2  determined  from  the  new  coeffi¬ 
cients  and  foF2  resulting  from  the  CCIR  coefficients.  Also  indicated  on  the 
figure  are  values  of  foF2  deduced  for  about  the  same  time  period  from  the  data 
collected  by  the  ISS-b  satellite  (Matuura,  1979).  While  differences  between 
observed  and  mapped  values  of  foF2  exist,  the  values  of  foF2  deduced  using  the 
new  coefficients  represent  the  observations  much  better  than  do  the  values 
determined  from  the  CCIR  coefficients. 

Further  validation  and  verification  of  the  values  of  foF2  determined  from 
the  new  coefficients  are  certainly  warranted.  This  is  particularly  desirable 
for  regions  of  the  globe  that  are  inaccessible  to  routine  ground-based  iono¬ 
spheric  soundings.  However,  it  does  appear  that  tne  new  coefficients  afford 
an  opportunity  to  determine  values  of  foF2  on  a  global  basis  that  are 
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consistent  with  observation  and  with  the  current  physical  understanding  of  the 
mechanisms  that  control  the  F2  region  ionization  distribution. 


6.  CONCLUSIONS 

A  new  set  of  numerical  coefficients  to  represent  the  median  behavior  of 
the  global  variation  of  foF2  has  been  developed  that  is  consistent  with  the 
current  understanding  of  the  physical  processes  that  govern  the  F2  region. 

The  coefficients  enable  median  values  of  foF2  to  be  determined  at  any  location 
on  the  globe  for  any  month  of  the  year  and  for'  any  phase  of  the  solar  cycle. 
The  coefficients  that  have  been  described  in  this  report  have  the  same  format 
as  those  developed  by  Jones  and  Gal  let  (1962)  which  form  the  basis  of  many 
ionospheric  propagation  prediction  methods. 

The  new  coefficients  have  been  determined  using  observations  of  foF2, 
values  of  foF2  determined  from  the  time-dependent  continuity  equation  for  ions 
and  electrons,  and  values  of  foF2  deduced  from  a  constant  local -time  analysis 
for  the  low-latitude  F2  region  distribution.  The  coefficients  have  been 
determined  for  each  month  of  two  different  years:  a  year  of  minimum  solar 
activity  (July  1975  to  June  1976)  and  a  year  of  relative  maximum  solar  activ¬ 
ity  (July  1978  to  June  1979).  The  time-dependent  continuity  equation  was  used 
to  provide  values  of  foF2  in  regions  of  the  middle  and  high  latitude  iono¬ 
sphere  where  observations  are  unavailable.  By  incorporating  realistic  models 
of  the  F2  region  neutral-air  winds  into  the  continuity  equations,  along  with 
a  representation  of  the  geomagnetic  field,  values  of  foF2  were  determined  that 
include  the  effect  of  both  the  magnetic  declination  and  magnetic  inclination 
control  on  the  F2  region.  No  direct  account  of  the  movement  of  ionization  by 
electrodynamic  (E  x  B)  drifts  was  included.  This  does  not  present  any  limi¬ 
tation  of  consequence  in  this  study  because  for  the  middle  latitude  F2  region, 
the  movement  of  ionization  by  neutral-air  winds  is  far  more  significant  than 
movement  by  (E  x  8)  drift  (Benkhe  and  Kohl,  1974).  At  the  low  latitudes  where 
the  (E  x  B)  drift  is  of  utmost  importance  in  the  F2  region,  the  continuity 
equation  was  not  used  to  determine  values  of  foF2  in  data-sparse  areas. 

Rather,  a  constant  local-time  analysis  was  used  to  determine  values  of  foF2  at 
low-latitude  regions  of  the  globe  that  are  inaccessible  to  routine  ionospheric 
sounding.  It  is  expected  that  low-latitude  values  of  foF2  determined  from  a 
continuity  equation  that  includes  realistic  E  x  8  drifts  in  the  equatorial 
region  would  lead  to  a  better  representation.  Such  a  study  should  be  pursued. 
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The  new  coefficients  were  used  to  determine  values  of  foF2  at  various 
locations  in  the  ionosphere  for  different  solar  activity  conditions.  It  was 
found  that  the  new  coefficients  yielded  values  of  foF2  during  solar  minimum 
conditions  that  were  on  average  better  than  the  values  of  foF2  determined  from 
the  CCIR  coefficients.  For  solar  maximum  conditions£  the  results  of  the 
analysis  carried  out  in  this  study  indicate  that,  if  adjustments  are  made  to 
take  into  account  the  fact  that  the  data  used  to  generate  the  solar  maximum 
coefficients  are  not  truly  representative  of  solar  maximum  conditions  for  all 
months,  an  improved  representation  of  foF2  will  emerge.  In  regions  of  the 
globe  that  are  inaccessible  to  routine  observation,  the  new  coefficients  yield 
values  of  foF2  that  differ  significantly  from  those  obtained  from  the  CCIR 
coefficients.  The  limited  validations  performed  in  this  study  for  those 
inaccessible  regions  indicate  that  the  new  coefficients  yield  much  improved 
values  of  foF2.  Further  validation,  using  data  such  as  obtained  by  the  ISS-b 
satellite,  is  needed  before  the  new  coefficients  can  be  recommended  as  a 
replacement  for  the  existing  CCIR  coefficients. 

It  is  necessary  that  further  work  be  done  to  improve  the  solar  cycle 
variation  of  foF2  determined  using  the  new  coefficients.  As  was  indicated  in 
Section  5,  the  values  of  the  new  coefficients  for  periods  surrounding  solar 
maximum  conditions  yield  values  of  foF2  that  are  not  an  obvious  improvement 
over  the  results  obtained  using  the  CCIR  coefficients.  The  new  coefficients 
must  be  adjusted  to  account  for  the  fact  that  for  some  of  the  months  (parti¬ 
cularly  July  through  December)  in  solar  maximum,  the  coefficients  were 
determined  from  data  that  may  not  be  truly  representative  of  the  maximum  solar 
activity  conditions.  The  relationship  chosen  to  represent  the  solar  cycle 
variation  of  foF2  obviously  will  influence  the  results  obtained  for  predicted 
foF2  values.  The  high  prediction  errors  for  both  sets  of  coefficients  for 
some  of  the  months  (see  Tables  7,  8,  and  9)  illustrate  the  fact  that  equation 
(8)  does  not  provide  realistic  predicted  values  of  foF2  for  all  solar  cycle 
conditions.  The  method  of  predicting  foF2  described  by  Liu  and  Smith  (1982) 
may  provide  a  useful  approach  to  improve  this  situation. 

The  most  promising  validation  and  assessment  of  the  new  coefficients  may 
be  their  use  in  existing  ionospheric  models  and  radio  propagation  simulation 
techniques.  It  has  been  mentioned  many  times  throughout  this  report  that  the 
new  coefficients  have  the  same  format  as  the  current  CCIR  coefficients.  Any 
ionospheric  models  or  radio  propagation  preo.ction  techniques  using  the  CCIR 
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coefficients  (CCIR,  1982b)  or  the  coefficients  of  Jones  et  al .  (1969)  can 
readily  be  adapted  to  using  the  new  coefficients.  Results  of  ionospheric 
modeling  efforts  and  radio  propagation  simulation  studies  obtained  using  both 
the  CCIR  coefficients  and  the  new  coefficients  would  prove  to  be  invaluable  in 
the  final  assessment  of  the  improvement  afforded  by  the  new  coefficients. 
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APPENDIX  A:  MERIDIONAL  AND  ZONAL  WIND  MODELS 

USED  IN  GENERATION  OF  NEW  COEFFICIENTS 

This  appendix  provides  the  meridional  and  zonal  components  of  the 
neutral -air  winds  that  were  used  in  the  solutions  to  the  time-dependent 
continuity  equation  for  ions  and  electrons  in  this  study.  Each  of  the  12 
figures  in  this  Appendix  (one  for  each  month)  shows  the  two  components  of  the 
neutral  wind  in  both  hemispheres  for  solar  maximum  and  solar  minimum 
conditions.  The  winds  were  derived  taking  into  account  values  of  foF2 
observed  during  the  specific  time  period  shown.  The  manner  in  which  this  was 
done  is  described  in  Section  3.4. 


neutral 


NORTHERN  HEMISPHERE  SOUTHERN  HEMISPHERE 


neutral  wind  models  for  June  1976  and  1979 


NORTHERN  HEMISPHERE  SOUTHERN  HEMISPHERE 


5  and  1978 


NORTHERN  HEMISPHERE  SOUTHERN  HEMISPHERE 


neutral  wind  models  for  November  1975  and  1978 


neutral  wind  models  for  December  1975  and  1978 


APPENDIX  B:  TABLES  OF  COEFFICIENTS  TO  REPRESENT  THE  GLOBAL 
VARIATION  OF  THE  MEDIAN  VALUE  OF  foF2  FOR  EACH 
MONTH  OF  THE  YEAR  FOR  TWO  SOLAR  EPOCHS 

This  appendix  provides  the  listing  of  the  numerical  coefficients  derived 
in  this  study.  Numerical  coefficients  for  the  months  July  1975  to  June  1976 
and  for  July  1978  to  June  1979  are  presented  in  tabular  form.  These  coeffi¬ 
cients  are  used  to  determine  the  value  of  foF2  anywhere  on  the  globe  for  the 
month  indicated. 

The  general  form  of  the  numerical  map  function,  &(A,  6,  T),  is  the 
Fourier  time  series: 

H 

fl(A,  e,  T)  =  art  (A,  6)  +  £  [a.  (A,  6)  cos  j'T  +  bi  (A,  6)  sin  jT]  (B.l) 
o  j=i  J  J 

where: 

X  =  geographic  latitude  (-90°  _<  A  j(  90°) 

9  =  east  geographic  longitude  (0°  0  _<  360°) 

(0  in  degrees  east  of  the  Greenwich  meridian) 

T  =  universal  time  expressed  as  an  angle  (-180°  <_  T  180°), 
and  H  =  the  maximum  number  of  harmonics  used  to  represent  the  diurnal 
variation. 

The  Fourier  coefficients,  a,  (A,  8)  and  b,  (A,  9)  vary  with  the  geographic 

J  J 

coordinates  and  are  represented  by  series  of  the  form:  -90°  £  A  <  90°. 

K 

a.  (A,  8)  =  ir  (^»  ®)»  j  3  1»  2,  ...»  H,  (8. 2a) 

J  k=0 


bj  (*.  0)  3  It  U2j-l,k  °k  9^’  J  a  l*  2 


i  •  •  •  >  H « 


(8.2b) 


The  particular  choice  of  the  functions,  Gk  (A,  8),  is  determined  by  specifying 
the  integers  k  (kg,  k^,  k^,  ....  k^,  ...,  kg ;  kg  =  K  where  i  is  the  order  in 
longitude.  Therefore,  a  numerical  map  can  be  written  more  explicitly  in  the 


Q'A,  6,  T)  «.'£  uo  k  Gv  (*.  9)  +  £  Ccos  1T  E  u2<  k  Gk  C) 
k=0  j=l  k*0 

*  sin  jt  £  Uy.,  „  t».  «)] 


0.3) 


Tables  B.l  to  B.24  present  the  coefficients  of  numerical  maps  of  the 
monthly  median  values  of  foF2.  These  coefficients  are  the  values  of  Us  k  (see 
equations  (B.2)  and  (B.3)  that  define  the  function  &(A,  0,  T)  of  the  numer¬ 
ical  map  for  the  indicated  month. 

The  month  and  the  sunspot  number  are  indicated  at  the  top  of  each 
table.  The  coefficients  are  listed  according  to  the  following  groupings: 

I.  Main  latitudinal  variation.  Mixed  latitudinal  and  longitudinal 
variations; 

II.  First  order  in  longitude; 

III.  Second  order  in  longitude;.  .  . 


IX.  Eighth  order  in  longitude. 


Table  B.l  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  January  1976 
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Table  B.l  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  January  1976  (Continued) 
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Table  B.2  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  February  1976 


HARMONIC 


.*2$25*m«oi 

.14)471201*00 

.101551001*0? 

-.*11171101*01 

-.*•175*501*0? 

-.170515101*0? 

.*424*415102 

.t«*msoi*oi 

-.17*11*201*02 

-.1700*9102*05 

.«4**44*9(«01 

*121110101*01 


-.121*70*01*00 
•2**5**501*00 
-.115101501*00 
—***2**tQt*01 
.20)177*01*0? 
.1215*1101*01 
-.7*71*5*0? *02 
.mis«*ot«o2 
. 100010002*01 
-.7111*1501*02 
-.«02*7|*0C*02 
.1***11501*02 


•«)*5**?t0f *00 
.501171201*00 
•11*20I10(*01 
-.2111**101*00 
-•**1411002*00 
-.1)5**0502*01 
.••1*11*01*00 
•11341*201*02 
-.51*50*201*01 
-.270*00702*01 
.)a«77*50l*01 
.***•**102*01 


10*1*1*02*00 

10**2*001*00 

711**1*02*00 

)***5*00(«0i 

71*1*1501*01 

1*7511*01*02 

.115101002*02 

125077*01*02 

*1**1*001*02 

2**755702*02 

21517*001*02 

105755201*02 


.7*51*2101-01 

.201)1*002*00 

•11150*001*01 

.7*12*2)01*01 

.21117*102*01 

•i*m«ioi*o2 

•*7)100001*02 

*507172201*02 

•*7*5*0701*02 

.*02*07502*0? 

.17*574)02*02 

.•1*0*5101*01 


'.50)021*01-01 

.72**1*102-01 

•125211*01*00 

>•1**22*002*01 

.*2*1*9201*01 

.50*1*1*51*01 

.mt*i*i>t*o? 

.4)7*7*151.00 
.2*7*0**02*02 
>.  **5100701*01 
.10*170)02  02 
.**l?2**04*0t 


-.***•01*02-01 

•simitotoi 

-.111175*0105 
-.5015)1101*05 
• 1*1012*02*01 
-.1)1711)01*51 
-.17)10**0201 
.525**2*02*01 
-.1)017210202 
.57*775)0101 

.m*o*?o«ot 

-.21*5)5*04*02 
. 151)7**01*02 
-.120555001*02 
-.7)5000*0102 
«7*l42%*0AO2 
-.**405*702  0? 
.15011101*02 
•0**l)050l*52 
-.072*1**01*02 
.21*00*1020? 
-.7510*000101 
-.)7*4»07020I 
.*270*7 200*01 


.224)25201-01 
.1**517201 *00 
.1*2)57101 *0v 
.1***55501*31 
.tumooioo 
.l*«*5010i*ci 

.nomtcioi 

.17*511104*01 
.tumioioi 
•  7*110101*01 
. ) 72  If l|Ot *0 1 
.1*17**994 *01 


•225111001-01 

.5811)001*01 

.702757401-01 

•1054i?00(*0* 

»124«4«5bi*91 

•  2M7jt*0<*3l 
.**7471001*03 
.7U7UT0I-C2 
•l*ft)157ul*02 

•  227*4*70  *02 
.5 1202**02*01 
•111*15*01*32 
.*3*2)4 >0k *02 
•12*51**01*0? 
.1*141270*  «02 
•*7*>,|501 *42 
.•IW«  5 101*02 
••*•14*1*01  *C  I 
.175*45001.42 
.13015**01  »C* 
.17714*101*3? 
♦*%*31»7C5*0* 
.53)521)01  01 
•481510531.0! 


b\  ' 


Table  B.2  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  cf  foF2  for  February  1976  (Continued) 
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Table  B.3  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foFb  for  March  1976 
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Table  B.3  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  March  1976  (Continued) 
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Table  B.4- 


Coefficients  to  Represent  the.  Global  Variation  of  the 
Median  Value  of  foF2  for  April  1976 
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Table  B.5  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  May  1976  (Continued) 
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Table  B.8  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  August  1975 


HARMONIC 


414)Miai«01 

•miittoio) 

.iwmnwi 

>.914070101*02 

,409010901*0? 

.1011*0001*02 


570111*01-01 

00t»»(<01 

10)71*101*00 

«l*31))Qt*01 

>}toi*e«of«e& 

mnotONi 

mm*oi*oi 

,10001*001*0) 

,notonoi*« 

,1079*7101*01 

771*17041*01 

,404)10401*02 


•22*1*1101*00 
•1**17*001*00 
-.111101*01*00 
-.5*2*10101*01 
-•1*5*7*101*01 
•5111*2 *tf*0t 
.10***1251*02 
•11**1*207*02 
-.*•1152301*01 
-.2 11707*01*02 
•21**47*04*01 
.117057002*01 


jrmiHHi 

>.|**2H»6t-V# 

■•imMlfl  *44 


2*1 *1*11-01 
.|K9H*Tf*4« 

.|41*Mi*9*) 


.1  l*o**OH-?l 


-  .a**aef-4i 


.*7*071101-01 
>.*9***2901*00 
.2)lt*»10f*00 
.1*4 2*0*0 f  *02 
>.*mil20C*01 
>.711222101*02 
.2)155**01*02 
. 1**7*0*01*01 
■•111*90*01*02 
>.191115201*02 
.1*00*7701*02 
•512*07 101*02 


-.91***2002-01 
.179*00701*00 
.101*77002*00 
.***1*0*01-01 
.**01*7*01*00 
•  2«*OMiOI*01 

-.um*iei«oi 

.?T\251S0(*01 

-.101**2201*02 

-.11)722*41*01 

.t*5on*CI*02 

-.7*0290101*01 
.1000*5*11*49 
»27*099**l  *01 

-•**077*701*40 

-.17*111907*03 

.timtt(v*0) 

.***79**07*41 

.*12299*07*47 

.11921*901*12 

•«50*0t*00t*4t 


-.*21199707*00 
— *«12I*I0(-41 
.*0*10*Mt«40 
.12*170707*01 
•  tllltlMMl 
-.109**7141*41 
•*lt**mi*ii 
.1*94*) *01*41 

»  7|  707*007*01 
-•UlioiTOCM 

-.m»**of«o2 

-.22447720* >4* 
-. *0*9*2  >44 -41 
.UimHI>92 
.27*217117*01 

»|922*U0t*O2 

-.im*)907*eo 


,*11077*01-02 
.141**1001*00 
•»1  *119*101*01 
•2*0***) 01 *01 
.0*10*7001*01 
.0209*0101*01 
•112**2901*02 
.(1)72)*0l*n 
.119722*01*0) 
,**00*0901*01 
.12212*407*01 
.*1*715)00 *01 


.27**21501*40 
.♦♦*«* 5*01-01 
•**1112001-01 
.* 1*21 5*00*00 
••93*15*0*1-01 
.21210)101 >00 
-.l«1l0*?4t >01 
>•0**771001*01 
.19*979101*01 
>,)**0)1I0(*01 
,4l9*)1*0t«0l 
.502517*01*02 
••IT2M070l«0t 
.10*12*100*02 
>•0*00*0*41*02 
>•  1027*0101  *02 
«1912)**0I*42 
>.1*2120000*02 
•11*211701*0) 
.991104400*42 
>. 11*7*7500 »0I 
.050*07900*01 
>.» 02*7 1*00*02 
•  .17*44*500*42 


>1 **1*1151-42 

-♦1*151 l»0!*4l 
,17*01 1*00*41 
>tMU»74*1l 
.**9 7* 7x70*41 
•*471*1401*41 
.175*^0100*41 
►I27ii5^00*«4 
-.155*77101*01 
.10*140500 •«> 
nmitmM? 
-. 11*11*3*1 »*J 

.nn**2o*«oi 

.ltt*9  K4I  »*} 


.lUttiiiHl 
-.itmitJi'M 
. 1212*70*1 ‘0* 
.5*7797941 **M 

.si 

-.44  7lt*MI  *24 
♦.TMtiviMl 

,iuttimm 

.4«*2**«4l-4l 
-.1 0159*102 *92 
—  14T««U14»44 
.imttitHi 
•.irntmis". 

.11*41*242  >-» 

•tunmiHt 
.123739  *01 *91 
.*.*7794*01  «44 

-.V>u  *r*tt-9i 

-. 141227101-41 

.?-V7*4Ut*42 

.ItUMilHl 

-.4»»5**7t1-tl 

.n»ftttM.K 

-.507*44401 -*t 

••uu*«u»-u 

-a«*2*55|l-4| 

-.2-45 14*11-41 

.*itioini-4i 

.4||574*t1-1| 

>.tt't»i%M  -t; 

.tt»9U4)t-7t 

•  -1 **  9  440 1*<  t 

>a*o*moi-oi 

-.110**1701*00 
.27*2*1001-0) 
.100404002*01 
•«12*00*t0l*01 
. 111*1*941*)! 
•2**1)2191*41 
1521***01*92 
>.21*0)5991*01 
, <*4105*04*02 
.•7*711001*04 
>•51004*191*01 


.10*21)007-41 
•.*041*1701-91 
••217950161*90 
-.20107*707*60 
.1*4101007-91 
.100174001*91 
.159120501*41 
-•775070702*01 
-•(1771)101*01 
.7744*0501*91 
,H*  11**01*01 
*♦204740711*01 


-•2«149»teo*tc 

-.*741*1101-01 

*4it»t«m*» 

» 1149044*0*60 
•  40|01>4l4t*04 
.I4«U-*U«*44 
.474*1 1*41 *40 

-.t*12*»101-04 

-,74t)S«><1*t4 

-.197(5*1(1-01 

■  '■ 

.1551194*1*04 
-.2514474*1 *09 
— ♦  1*  *'*54  ioi  *04 
-.2*721 0701*44 
-.**»*< till -44 
ai;t*t(»c*»l 
•  UttlitllMl 
.741970111*44 
.1*2*71104  5)1 

..mimiH) 

.*72010191-42 

,»»«7*3*tl-Ot 

.ijmntt'it 

-a«49*:»it*tt 

♦.imttNi«'i 

.110*11*27-41 

-.11*-  >>WH 

-...tnjtttm 

.10)005002*40 
-.500*14541-01 
•)IS*4|)tl*04 
atmM(*M 
-.12«CVU0I*9L 
.140151714*01 
-.1149*4  lit *94 
-,*4*?4049t«^4 
.00*044511*41 

«.iiiimoi*oi 

•14444*092*01 

•14)747491*40 

—i*iiuio2*t: 

.1)1404*92*42 
-•121)9*091 *42 
.mu*tit*M 
.10x4*1192*41 


-.*07210144 *40 
•  Mt)t4Mt**t 

vft*4U+t}««9 

stmumui 

.102*51)7*4) 

-.M47*«*2**0 

**»)*«  «*«*« 


•201141241-41 
.102124100-41 
-.4*1107201-01 
,)f 1)20041*9* 

—  )5)fim4*9ft 
•♦102117791*90 
.11ft0»4«0t-tl 
••4*1411002*41 
♦1090*1402*01 

♦inm*«(*«i 

-.141*4  Id  4 -*i 
♦123044*01*42 
•♦11244-1)41*01 
-.*t«m*e 2*oi 

104**4404 *02 
•119241267*41 


.1*2115104*4) 

•  )i)t(:M;«o) 

.417)45*94-01 

atinjiso-ot 

.074*5*154*90 

•  10*t»01€4*«4 

.41(522411-01 


.m»47tt0-«1 

,4tM**l«>t 

-.♦r\t:*44*-«l 

,799721751-0 

.t»*»*32*i-o; 

.m*i«*oo-ot 

-•42«*Q\t«t-42 

-.51*11)1*2*41 

.13 *144791 -41 

•.1419*5001*04 

-.1*790*140-14 

.t«57»k40| -ft 

-.<175*0121-01 

.101571141-11 

-.lrt1i*12»-». 

-.7*471*105-19 

.14*5*9501*4 

.tl)5**|7i-  ; 

sittUtttiMl 

-.MKHUl-f 

sttuum-os 

5.)  *v*2  ?♦**-** 

-.5*1*90101  St 

.«4«Vt?7*t.?J 

V.  j>‘W 


Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  August  1975  (Continued) 


Table  B.9  Coefficients  to  Represent  the  Global  Variation  of  the 
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Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  September  1975  (Continued) 
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Table  B.10  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  October  1975 
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Median  Value  of  foF2  for  November  1975 
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Table  B.ll  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  November  1975  (Continued) 


Table  B.12  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  December  1975 


Table  B.12  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  December  1975  (Continued) 


Table  B.13  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  January  1979 
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Table  B.13  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  January  1979  (Continued) 
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Table  B. 15  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  March  1979 
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Table  B.15  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  March  1979  (Continued) 


.tmiiHHi 

.  7*012101  *oi 
••ittnutroi 

-  niimtHt 

•.tMUUHl 

.MMSK9K( 

..OtUlHHf 

•.U’MMH) 

.MwitcHt 

.**•  *l<Ut*09 

.lKI1IU44i 


4i«ui4i«m 

nlMIHUMl 

.looutv-oi 

•.(IllilUHt 

-.om»*w*i 

.mimiHi 

•.iHummt 


-.MOIMI'OI 

tit*u«ooi*et 

.119471111*12 

-.MIOUMI'M 

-.mostiaH) 

•  Hl*MOf*47 

•.tnnuoioi 

.HITUUhSt 

<.um«48i<oi 

-.witua«|.9) 

.»K3tH«l«92 


•.nutiid-at 

•nimu(*oi 

.uioom-o 

.4IIIIU4M91 

••w**»*«9«oi 

.UHtMIMl 

.nnMUHj 

.llltHHHJ 

..UHMKIiU 


.itimtsKi 

-.«7t9*«l9}*C3 

•MttUKhW 

-♦to7k*to»*oi 

.tiettmNi 


77  (19101-41 
.(MmitHi 

•tUiMOHI 


HetauiHi 

uunwKi 

PJ»1|49l*4l 
.miMStt'Ot 
-.(*(*•  *«•« 
•.UJ4  *11^ 

.Ifu'fMHl 

.huiomh: 

.iTniTt«i4t 


-.mcioKi-ji 

•  110* t*l Jl *54 
.moioKhca 
>.4icta*)oi»oi 

,l*K10C»*Sl 

.V»>4*M9I«41 

>»me«oioi-o2 

'.i«n»noi*9i 

.itKawci>oi 

•iiM2nHi*oa 


.iocmt(Ki 
.loiaeiuMi 
tmU«V9|*09 
>•  1441**101-01 
.IKOOIOI-Ql 
•.1(627*409*61 
.44100U5f*01 
•I163?)l4t*ai 
•mmioi-oo 
•uno^oi*4i 
.tU^roio: 
.iHduaioi 

•  27*44091-02 
-‘•OMOIOI-Ol 

•  J-.*2O79T0*-O2 
.UIIH8CI  *02 


♦mmtai-oi 

*•01219(101 *C9 

.imuisNi 

-.n>M(oo4»oi 

4.1912*7101*02 

•1MM540I  *01 

-.(tieietooot 

-.•uiioweot 

.)»U*»UC»40t 


.72)0)401-02 
,IM2M«0I4I 
s 11 1714101*00 
.((1471969 *«« 
am()99t*oi 
*ni«t«toi*oi 

•.9*2*7**01*01 
•110172(01*92 
•  19HQ4*0t*O2 
'•12*1(0101  *01 
••*16*21101*21 
•**070*01*91 


t«M0210l*01 

2076*1101*06 


atiMlitHi 

-|*4(tl%*0*«0 

•.roHJU|-« 

sim*it«m 

.Itlituiu* 

•  (  *  1  **♦*  *  t  ^ 

••9**21**69-61 

••uj:*2M9«>e 

,*»»:***«>« 

•  n*7*2**«»C7 
•  .iiltUKVM 

.in«Hu»*u 

•ttmott-o 

•Ml*(t*t*««l 

•atutmwi 

•aiM-.iHf**,* 

.i*t>9»*7t«24 

•a»ft»mt*n 

•  s*:r:  Mnt-« 

.:***(»((( »u 

•t0i|**«t|>49 

..itimiTt.ii 

.7-««:u*t-ti 

.MiHIMHl 

•.mmtiH. 

—  ttOO*17*H 

.ut»4«nf*67 


.4**117(41-21 

.2UtMU2«61 

«tt****nt*n 

•aitutuM; 

aotuoftt 


.UU}i*iW} 

,{!»»•»«  -it 


•  ttdflMt-tl 


•  It  H*  "**!•*( 
.;»niu)i.cc 


.*»**»*«  1***1 

-o 

•9*»«*;;i<«4l 

•  .?U*»«*t**-f; 

Table  B.16  Coefficients  to  Represent  the  Global 
Median  Value  of  foF2  for  April  1979 
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Table  B.16  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  April  1979  (Continued) 
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Table  B.17  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  May  1979 
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Table  B.17  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  May  1979  (Continued) 
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Table  B.18 


Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  June  1979 
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Table  B.18 


HARMONIC 


Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  June  1979  (Continued) 
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Table  B.19  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  July  1978 
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Table  B.19  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  July  1978  (Continued) 


Table  B.20  Coefficients  to  Represent  the  Global 
Median  Value  of  foF2  for  August  1978 
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Coefficients  to  Represent  the  Global  Variation  of  the 
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Table  B.21  Coefficients  to  Represent  the  Global  Variation  of  the 
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Table  B.Z1  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  September  1978  (Continued) 
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Table  B.22  Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  October  1978 


Coefficients  to  Represent  the  Global  Variation  of  the 
Median  Value  of  foF2  for  October  1978  (Continued) 
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10 


11 


aiMimi-ti 

MinuitHi 

.UlUtlOH) 

•.ui*om«*4i 


,Mam90Hl 

•.*«s*nooi*t9 

.217*74)07.01 

.tmmoi'81 
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.I47||449f-0I 

-.191)42791*01 
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,|M«m9l*91 
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.14*124101*91 
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■.»7MUH*9» 
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-.)4tM449lf  1 
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■*9927*4991*91 
.1147*4499*91 
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•.949*21791*92 
.111429799*91 
••11*499*91*91 
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•971411091-91 
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.744474991-01 
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-.142947191*02 
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•*149%9I0| *92 
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•.*974*2191*92 
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2971*4401*92 
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.2HC4|)9l*3l 
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Table  B.23 


Coefficients  tc  Represent  the  Global  Variation  of  the  . 
Median  Value  of  foF2  for  November  1978 


HARMONIC  o 

1 

2 

3 

\S 

K\ 

0 

1 

2 

o 

4 

5 

6 

0 

.101527406*02 

-.576419706-01 

-.14722770E-01 

-« 57952410E-01 

.199157406-01 

-.678355006-01 

.196026006-01 

\ 

•  46155140?  *00 

-.24»«3int*oi 

.666241302*00 

-.413801806*00 

,158849006.00 
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-.590470206*00 

2 

•3576*090? *02 
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-.211793 10?  *01 
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i 
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.920373)05*01 

A 
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•62,02J506*01 
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-.143077206*02 
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I 

5 

-.2927*1406*01 

-.24304Q20£*03 

-.19406180E402 

-.307316506*02 

-.273880206*02 
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-.44)030001*02 

0 

•325925*0? 40 J 
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-.1*4394 10E >02 
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7 

••5055540? *03 
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.314421306 *02 
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• 

‘••'71333901*03 

.49782530? *02 

•10663I40E«02 

-*519377506*02 

.304701606*02 

-.208319006*02 
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9 

-.4  7041406*0) 

-,51455130?*03 

-. 565479905*02 

-.438486406*02 

-.377544006*02 
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-. 34841310? *61 

.228102106*02 

-.111237006*02 
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-.9113)8206*01 

u 

•401 40480? *03 

*1799-.i30£*0S 

,116216905*02 

.121880906*02 

.782767506*01 

-.192282106*02 

•247327205*02 

12 

-.118937806*00 

.151*2320C*01 

•19263030?  *01 

.132927206*00 

,940606906-01 

.26)1749^6-01 

.394800901-01 

f  1 

—  5833H50E*00 

-.224929406*01 

.13115490? *01 

—  234219206*00 

.223901706*00 

-, 714402406-01 

-.810060305-01 

Is 

-•114132406*00 

.!4988260E*00 

•30919140E*01 

-.26328150£*00 

.171728506*01 

-.48)640306*00 

.171927006*00 

1* 

-.127114906*01 

.  3900*1 80F*00 

,72366010F*00 

-.137304806*01 

-.913382806*00 

-.981923806-01 

.733028805*90 

16 

.162709206*02 

.184067001*02 

-•13586370E  *02 

-.350702006*01 

.367312106*01 

-.727449106*0) 

.591973306*00 

17 

•  25791630?  *01 

•23852520E*02 

.146632001*02 

-.107135706*01 

-.344044806*01 

-.448763)05-01 

.133203306*01 

18 

•2I956610?*02 

.343745708*02 

-.46310170?*02 

-.194730306*02 

-.112983306*02 

,103310706*01 

-.332998305*01 

19 

•60269440?*02 

-.*5454480E*01 

.258962206*02 

.943927506*01 

.168752006*01 

.8044)7806*01 

-.134)03406*02 

?n 

-•74991790E*02 

-.215240506  *02 

•  U929380E*03 

.144670306*02 

-.221424306*02 

.847030)05*01 

-.40900)201*01 

?1 

.444971706*02 

-.197289105*01 

.378476306*02 

.147110006*02 

.0)0045306*01 

.361100005*01 

-.236708)05*02 

22 

-.122077306*03 

36326530E  *03 

. *1069100?*03 

•19669360? *03 

.500945306*02 

,16927)106*02 

.316460705*02 

21 

“»316“'0540?*03 

-.U572640E*03 

-.177427306*03 

-.131467801*02 

-.566120006*01 

-.573303606*02 

.706897801*02 

11 

26 

•14904190?*OJ 

-» 17081 780? *63 

-.376123406*03 

-•  243948706*02 

.443663306*02 

-.239*79306*02 

.712173005*01 

21 

-.2211)2906*03 

.52976950?*03 

-.412232006*03 
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-.231329306*02 

,813131105*02 

26 

•  \4725760?  *03 

,833t6100?*C3 

-.126177906*04 

-.423046306*01 
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-.688680406*02 
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27 

»61009I30E*0  1 

.  50072070?  *03 

.422714706*01 

.699858706*02 

.531477406*02 

.157929105*03 

-.189269105*01 

2* 

••2)9141608*0) 

•  373*9J9Qt  *03 

.468082706*03 

.181041006*02 

-•411476606*02 

. 239192006*02 

-,1)3739405*01 

29 

.1270109* . *03 

-.50697410?  *03 

♦  75178740?  *03 

.480071506*02 

-.173672706*02 

.315443206*02 

-.10974I705-03 

10 

•95205160?*01 

-.800*441 '!*03 

. 1653337CE *04 

•47681960E*03 

.214651406*03 

.933003406*02 

,047747705*02 

n 

-.490711706*01 

-.70295030E*03 

.423642606*03 

-.766391 806*02 

-.103231206*03 

-.168763705*0) 

.200097701*03 

n 
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-.20501440?*05 

-.213380706*03 

-.546383406*01 

.1)6979506*02 

-.078119706*01 

-.624726105-01 

31 

-♦H042120l*0j 

.262463605*93 

.403592806*03 

-.206457606*02 

.763306706*01 

-.138605606*02 

.122188206*02 

3* 

-.70291)906*02 

.2?437?10E*0) 

-»73326*10E*03 

-.  192860306*03 

-.105373406*03 

-.431239106*02 

-.203107605*02 

31 

♦  14844550E*<}1 

.32223160?*01 

•15397020E*03 

.316478506*02 

,374707006*02 

.821954605*02 

-.777116105*02 

3* 

.24O63610?*00 

-.712 192  40?  •-.•4 

.87t6?020f-0l 

-780970306*00 

,756127006*00 

-.689520705-01 

-.15)141805*00 

17 

.150514106*00 

-.142*6540?  *0t 

-.41241870E*00 

-»  7 J2977806  *00 

-.741712806*00 

-.573485)06-03 

-.276247405-01 

31 

•1 3132440E*0l 

.721702601*00 

.101224406*01 

-.213914106*01 

.118361806*01 

.809146906*00 

-.153414705*00 

19 

.606531206*00 

-.12724820?  *00 

.241*62906*01 

-.643310006*00 

-.196176306*01 

.44)527906*00 

.1)506)405*01 

*0 

-.164»J630E*G1 

« *02471*06*02 

-.9Bf485?0E-0t 

-.124341106*02 

.311741006*01 

.372490906*00 

.414606105*00 

*1 

-.1089*1806*01 

•  7f 703500C  *01 

•11981680? *02 

-.444009506*01 

-.125325206*02 

.2)4954006*01 

.61)404905*00 

12 

-•S*' 486?0?*01 

•«11412330?*01 

-••14753506*02 

-.116781006*00 

-.783289706*01 

-.529744006*01 

-.112762405*01 

A3 

-•315104901*01 

.205072101*01 

-.2473*6906*02 

.3il53.'l06*01 

-.254901006*01 

-.1)6227006*01 

-.9919)1301*01 

64 

.21122180? *02 

■.209113706*02 

-.304791106*02 

.706522906*02 

-.210040306*02 

-.)097*9'06*01 

.110963705*01 

ill 

*1 

.2485*4601*02 

-*«  1684  *050?  *02 

-.544120001*02 

»12>'2160f *02 

.085932606*02 

-.15)5)9101*02 

->10!66890I*Q2 

*6 

•  UT04llQ5*07 

-.730964201*00 

.245*22306*02 

.123368006*02 

.8)6244106*01 

.117128906*02 

.777217101*01 

47 

.111170706*01 

.2976 J610? *01 

. 7011 JOOOt  *02 

-.•16173006*01 

. 2)045)206 *02 

.22 3014)01*01 

.200901106*02 

49 

-.26449440?*02 

.478215201*01 

.163(9890?  *01 

-.123966906*03 

.215178506*02 

.1)9549806*02 

-.612990105*01 

49 

-•406894)0**0? 

.398941801*01 

.12708290? *03 

-.216217606*00 

-.118)05206*01 

.312816106*02 

.298441005*02 

10 

-.702742204*01 

.284743>OE*0l 

-.18607040? *02 

-.673808 1 0?  *01 

.047949006*01 

-.872136306*01 

-.681410405*01 

n 

-.112151806*01 

-•100024906*02 

-.56142080?*0 2 

-.220230506*01 

-.190427906*02 

-.187188106*01 

-.111112105*02 

52 

•2*183450? *01 

.105601606*02 

-.133416402*03 

,72 1 2  *  0005*02 

-.422489706*00 

-.132514306*02 

, 7?7»9UOI*01 

31 

.406278101*02 

,111402101*02 

-.86581660? *02 

— mcmo?*02 

,6895)7701*02 

-.19761)706*02 

-.227642605*0) 

94 

-.21819590? *00 

1 3676140?  *00 

-. 223054906*00 

,  U0701406  *00 

.14940)006*00 

-.316140006*00 

-.901)11)01*00 

51 

.597748206-02 

.594022301-02 

.129834006*00 

-.112990606*00 

.169008706*00 

.992642405*00 

-.211140601*00 

54 

•  16108400?  *00 

.117071806*00 

-.13505360? *00 

1 26691 510! *00 

.11444)406*00 

-.1)992S)oe*O0 

-.445907101*00 

17 

.6183*010?  *00 

-.2*722500?  *00 

—  »  31959010?  *0" 

-. \9  4*01506 *0C 

-.7540*0206-01 

.  1 1985)406  *01 

-.174)4)701*00 

IV 

99 

»?60I3840?*01 

•»3508J?00?*00 

.tU36190?*00 

-.177018206*01 

.240142106*00 

.172444201*01 

.1)1)47401*00 

54 

-.741*8230?  *00 

.19U?U0?*01 

-.17929030?  *01 

-.644423506*00 

-,24414 )|06*Ci 

-.111171805*01 

•11161)401*00 

40 

-.23422200?  *01 

-.414265001*00 

-.724007706-91 

-.167741906*00 

-.294007606*00 

-.48302>2O!*OO 

.u?m40i*oi 

61 

-.2t7\l  840?  *01 

. 14364300?  *01 

.13227430? *01 

.2919370O?*OO 

.471406106*00 

-.246)71801*01 

-.724179401*00 

62 

-.13|72540?*01 

.102812806*01 

,101564801*00 

.231282806*01 

-.1244  7  4*Q?  *01 

-.)063«)70<*01 

-.611)46)01-41 

43 

-.387303701*00 

-.242287801*01 

.651131801*01 

.180445801*01 

.440990106*01 

,413090901*00 

-.1101)4)05*01 

64 

-*4t»m*0?“02 

.101153201*00 

.26T9947Q?-01 

.43 160040f-02 

.24517420<-9l 

-.11)3**106-01 

.1719*4101-01 

65 

-.117907106-01 

.4274441CI-01 

.100182301*00 

-.11112*106*00 

-•445433*06-01 

-.110944006-01 

-.7»)71M0|-O2 

V 

44 

.717737101-02 

.204731201*00 

-.569199001-01 

.5*7064606-01 

.10)129206*00 

-.301001006  *01 

-.411111605-01 

6* 

.100227701-01 

.123119201*00 

.279614801*00 

-.181154506*00 

.  W4>M90?»00 

-.22*112405-01 

.136714005*00 

VI 

69 

-t93394*4C?-01 

.121373401-01 

,784664101-01 

-.175177906-01 

.*430)100101 

-.2)4**0008-01 

*14*2x01-01 

49 

-. >11*41201-01 

-.277764008-01 

.116832801*00 

-.229474101-01 

-.142720106-01 

.4*71*1205-02 

« .62104) lOl -0 l 

VII 

70 

.100420*01-01 

.234194301-01 

.661594201-01 

-.148*46501-01 

-.m?onoi-o; 

*. )I04071Q|-01 

.47M146CC-0I 

n 

-.7014441CE-01 

,467040401-01 

.311461 10f-01 

-.417174006-02 

•  l44l *  2701  0 1 

-» |**2 1)808-01 

-.49210)401*01 

VIII 

72 

-.611101401-01 

.123360408-01 

-.473103701-02 

-.14629)406-01 

•.238440101-0) 

- .20Q40X3I-32 

-.124113*05-01 

71 

-.129)42106-01 

,130032901-01 

.224190301-01 

-.101017401-01 

,11 4*49061-) 1 

-.1X0  21  >08-01 

127*01201-01 

IX 

74 

-.103748001-01 

,299*4410?-01 

,273697201-01 

-.2»»092I06-01 

-.4*4910001-02 

-.12**91905-31 

.17061)601*41 

75 

-.374734001-01 

.419117101-01 

.463647701*01 

.4O»*510l  *01 

.IJ'O.'.MOI-Ol 

.»4»00M08-«t 

.... 

112 


-ft.. 


A 


* 


.to*962?o(-oo 

.4M(«»9I«40 

>.5000*419100 

>.926*9460(01 

uiosmroioi 

.17***1)0202 

.UH»)O(«0l 

•♦710T4OikO2 

•ioouosoio* 


.{•mtMHt 
-.156**9*0(00 
-♦  712420001  00 
.itmuocoi 
.201710901*01 
-.102095201*01 
-.11011*401*01 
•10*215*01*01 
♦17100*701*01 
-.2105*2101*01 
-.111501001-01 
•*5*1*5*01*00 


-.70642*60  (  02 
-.00*090200-01 
-.**•0*0*02-01 
.111150101*01 
.1616*120(01 
•.*10*0*101*01 
-.*10010001*01 
.9170*2*01*01 
.7192*5901*01 
-.1*01*0100*01 
-.219170500*01 
-.*1*77*001*00 


-.27*7*7*00-01 

-.111127101-01 

.9179*1201*00 

•971*1*001*00 

-.0090*2200*00 

-.;i2070*0f*0l 

-.202*1*500*01 

.*10*7**01*01 

.711*25700*01 

-.17101*501*01 

-.*11702*00*01 

-•**2010*01*00 


111220501*00 
1111**700-01 
11010*100*00 
*9**0**QI *00 
117121*00 *00 
10*1*1201 ‘00 
« 1 10*01 of *«0 
74441(069  00 

i*mii*t«oi 

212711100*00 


.57*273*91-0) 

.101799101*00 

.*9*711*01*00 

.205*17000*00 

.12*01*001*00 

.iM*oo«ei*oo 

•.*0**09*01  *00 
.*29710101*00 
>.112*02201*00 
•I792064CI-6* 


vii  R 


.2*474730(09 

♦2*05)44(00 

.1141405(1-00 

.15215110(04 

-.2*1444731 *Q« 
.2104*1)01-04 
-.475047)4(00 
.ItllllOOf-OO 

-.13901)14(03 

-.14*5*424(00 

.1**44140(01 

-.11-7*1*51. *6 

•♦2*441411(01 

-»ll*7?l*0(-01 

•,7M»*»|4IO* 

-.(4  719*29(01 

♦27«9>mroi 

-•)IO*)l40iO| 

.1)41*4*9(01 

-.7*15*1141*41 

1**07*700-01 
122*0*701*00 
l ••♦*1500*01 
11721*901*01 
1510/5*01*01 
1)1292201*01 
,;071H10I*02 
5MI9MOO*Ot 
101*11101*01 
11000*001 *01 
1*5*12001*02 
910)71100*01 
*1*911901*01 
21790710(01 
1)0911701*01 
.*17129101*01 


.977711)01-01 
.79*1)4201-01 
•19*011102-01 
.1)11*1002*00 
. *1**17T0| *09 
•111*01101*0* 
.297921202-01 
.111*21*92*00 
.19*977501-90 
.19*00*191*99 


-.1*7111102-01 
.1  1*19102-02 
.|91*5*tC2-Ol 
.9*2l*»»CI-Ol 


.11*1**502-00 

-.7416***3100 


•.917219707-01 
.717910101-01 
■•17*172202*00 
>.21*091901*00 
•*«90:i902*00 
•.117*71202*01 
.11*971701*01 
.272*71501*01 
-.1*0171702*01 
.*99091102*01 
-.1)2*09*01*01 
-. 71*99*101-01 
*.  14*41120(00 
-*10*93010l-02 
-.9*1190001-01 
.171115)02*01 
•120172*01*01 
.91)«)*«C2*01 


.291.191102-01 
-.211*29*02-02 
•••*9^7)0*01-01 
.17*117)02-01 
-.2*1)71202-00 
•  11>1I*1*C2*00 
.20)971*02*01 
-• 199200902-01 
.l**"9**02*00 
•.27*1**202*00 


.2*0)9*202-91 

•.*41*71*04-11 

•♦***10*201-92 

.1**7*1201-01 


. 29**1*105  *30 
.1*00*1902-01 


.1570**106-01 

.99)199901-01 

.17*010901-00 

,1*4*4460200 

.9)127*201-01 

.999277101-01 

■.2921-5202-01 

.107171001*00 

.919999102*01 

•*193*1999*00 

•.1*9001106*91 

-.991119*01*00 


-.295977*09-01 

-•*1*21*101-01 

-•99M2019I-01 

•22)0>910I*0Q 

.999921101*00 

.1**111502-91 

.115*1*101-01 

♦.*«omioi*oi 

>.919519901*01 
>.729*41501*01 
>•*92217901*01 
.272*1*901*02 
.2715*2901*02 
.19717)001*02 
.21*929002*01 
>.799*19502*07 
-.195994101*02 
-.1911*4701*07 
.4*4147402*01 
.7*9110101*02 
.177291201*02 
.90474*101  01 
..47175*401*0: 
-.119**7102*02 


-.*14191901-01 
.727*1*401-02 
.1*0*7100100 
.122014202-01 
.5295**70201 
.741*024101 
•.*240*1792*03 
-.109  720*0101 
.279150102*01 
-.147224101*01 
.101*91402*01 
.2944*240101 
.**91*940101 
.19211 97CI -07 
>,4171*000103 
2*4*5*40201 
25*12440201 
.9449*4492 -Cl 


,14791170601 
-.2)9142102*09 
-•*l«»2990IO3 
•*2117*20(01 
.204119002*01 
-.21 71 74131*02 
-.076*2702*01 
•  *4040*4K»w7 
. 5*599693?*  31 
-.♦l*l*?t  3202 
-.12159040201 
.1*911940102 


.427471101O1 
-•10*92*102-02 
-•*11701201-00 
-.20924910203 
-.1245**101-01 
-• 1421)2901*39 
.4*147)402-01 

♦  )77*0770l*-*l 
.771*1*10101 
.221611201-00 

-.11177*43(02 
-.194507502*5. 
••20*27)20**02 
.20711)102  01 
.772097442-05 
•*11*21232*32 

.2974***4.2-42 

-.*4)01610201 
-•617**4401-02 
-.17412*00(01 
-.121**7101*32 
.25114)402- 31 

♦  122*17442  02 

♦  12*12102-02 


-.*60410102-31 
.1022)0)1-31 
-.10**9140103 
.1)1*71)92*03 
.I41949l0l*)i 
••11*1*7402*03 
.*41245*01*09 
—.160967901-01 
-.10*90113101 
•17740**31*01 
-.14)94169101 
,92*mtoioi 
♦111*1210? *02 
-,)))9«2412-vl 
.132072502-01 
*11021 40t*U 
•.11144*40(02 
♦  2K*HO1*01 


1224*72)1-01 
15141413(01 
17*91650*01 
4*44*1*3(01 
211)1*191 *03 
*1*11110101 
4027*990(01 
11)5*1*11*03 
4121)07?(*:9 
77*79 1*3( >00 


ummioi 

1*4615*0101 

1(1)7573^-03 

1191)479101 


.2144)120102 
.1044*6334-32 
»  5*053*05— Jl 

•  •150*7)I9(*)3 

♦  *477577.(01 
■.H2*l**wt-*1 
>.(4l7**10(-w3 
.1*1*1371103 
■•U*27*4‘(-21 
•6)l)*760(-*l 


. *. l))*))l-32 
.viii*imo5 
»MQ4«}0(-:i 

♦  i:t(«n<j<-n 


•  *(7***035Ol 

-•<2*i**4CI-01 

♦73*7*190(01 

,)7*7|55«(-n 

♦  150111)5(01 

->244*1 3(3101 

.771*2*64(01 

.10370)10(01 

.l»l*?|»9(01 

-.155(5**110) 

-.*4-61*12(05 

|*7t*»*0(-i; 

Table  B.24  Coefficients  to  Represent  the  Global  Variation  of 
Median  Value  of  foF2  for  December  1978  (Continued) 


HARMONIC 

4 

5 

6 

H 

7 

8 

9 

10 

11 

12 

0 

.141*21901-01 

•**4010*01-01 

•27101*201-01 

*.200)4*«0(-01 

-.*10*1)231-01 

I 

■K  i  Jttrl 

-.201400201*00 

****40*01  *00 

-.70*07*201-01 

.2*179*401*00 

.2071*110(499 

2 

-.201441201*00 

-.* *791**0f *00 

•«714*4*90C*00 

.*2*117401*09 

-.117*44191*53 

1 

»MHJ090f*flC 

.9I1710*0|«01 

.11**44*0(*01 

17*4*7*91*01 

-.21442*10(401 

* 

■trisEi 

.*44171*01*00 

•2*19)0*0(401 

.*mrvio(»oi 

,.2**141)0(401 

.234*7929(401 

I 

1 

»*T)*3340f *32 

.292490104*01 

-.212*79701-02 

-•*31*1*901*01 

.17207710(401 

•10*12**0(492 

1 

ft 

-.904111401*01 

.747*24400*00 

-.9*010410(401 

-,142*2*401*02 

.9*0***10(401 

-.7444*700(01 

T 

-.4iumoi*ot 

-.199117001*02 

.*09*1190(402 

•2*14*1*01*01 

.*4947500(401 

-.20914420(431 

« 

•*11041*0(441 

-.1**174400*01 

.5*949920(401 

.21*440*01*02 

*.*17*7400(401 

.4219)940(491 

1 

TTTTTTtTTl 

.21411*401*02 

-.112**700(402 

.)244l*00(*0t 

,.47977070(401 

.17140990102 

10 

.199*10401*01 

-.20974740(401 

-,101179*0(*0l 

. *4341 *706*00 

-.1**»4500(401 

11 

-.492771901*01 

.171 410 *0(401 

-.117102 902*01 

.*9491*001*01 

-.94I2*«70(*01 

u 

-.2*1*24001-01 

•9*2A7*90(-Ot 

-*)U7U*0(-02 

-.711*2*001-02 

-.1)7)** 732-01 

11 

•400*72*01-01 

•11)13*701-01 

-»114*1*50(-01 

-.*4*)*tl0(-31 

14 

.110549701*01 

.4427*4901-01 

-.129910*0(401 

-.1*2*51 *0f *00 

, *9*077*01.00 

•14*10*90140) 

11 

.4*19TT10I*00 

.499144001*00 

-.1*77*1)0(400 

-.*7441*101*00 

.17710*10(400 

,14419040(403 

u 

Iff  TTTTl 

.107744701*01 

.40*71)101400 

.9)40***0(400 

-•2*40)lt0(>00 

-.7*1*10101*09 

IT 

-lift  2't'i 

. 412419200*01 

-.14994*10(401 

•**9072*02*00 

-.1*0441501-01 

-.1279*3201*00 

la 

Ef . it  .It  la 

-.212471200*01 

.20912740(407 

.1*9272*01*01 

-. 3*942*101-01 

-.104201*32*11 

1ft 

*.T4ft4fttlO|*01 

-.*14*79*00*01 

.91104110(401 

.524*19902*01 

-.121*9**0(401 

•,9944*4)01*01 

« 

—.21474*701 *01 

-.40*51*001*01 

-.11*79040(402 

•.921040401*01 

.1*1**090(401 

.mi)*)0f*0l 

tl 

-.114714501*01 

-.11*177901*02 

.119*1100(402 

-.122)97102*01 

.1*109*40(401 

,)914*M<1?*91 

» 

.411 144108*01 

.109*74100*02 

-.111*11*0(401 

-• 17I49(*0(*02 

.2**7*9001*01 

-.55149*901*01 

11 

li 

.444114001*01 

.11**79*00*02 

-.294**0*0(402 

-.21050*40(402 

.42700470(4*1 

•!41499Q0(*Q? 

14 

-.150944401*00 

•2*9*19701*02 

.479*79*0(402 

.14217**0(402 

-.20014)10(401 

-.24*121001*92 

11 

.149441401*02 

.10*10*100*01 

-.*0*17270(402 

»*2l)l4*0l*0t 

-,*41*4*401*01 

11*41030(0! 

1* 

-•1114T010I *01 

-.11944*900*02 

.2**797*02401 

.*0477*00(402 

-•47***9101-02 

•it**4**)t*92 

IT 

-.110414401 *«7 

-•11*1*1*01*07 

.9**107*0(402 

.**>921*0(402 

-.14970713(402 

-.*4170**01*02 

11 

.TM71110I»01 

-.119*02*00*^2 

-,«19**T*0|4C2 

-.270147*01*02 

-•*21141701-00 

.*12**9231*32 

14 

-.114101401*01 

.92**2510(402 
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APPENDIX  C:  CONTOUR  TABLES  OF  THE  GLOBAL  REPRESENTATION  OF  THE  MEDIAN 
VALUE  OF  foF2  SELECTED  MONTHS  AND  UNIVERSAL  TIMES 

This  appendix  provides  examples  of  global  maps  of  foF2  determined  from 
the  numerical  coefficients  derived  in  this  study.  Shown  are  contours  of  foF2 
in  MHz  for  0600  and  1800  hrs  universal  time  for  the  months  of  January,  April, 
July,  and  October  for  solar  minimum  (1975/1976)  and  solar  maximum  (1978/1979) 
conditions.  The  maps  illustrate  the  global  variations  of  the  median  behavior 
of  foF2  as  a  function  of  season  and  solar  cycle.  It  is  worthwhile  noting  that 
the  gradients  in  foF2  are  less  steep  on  a  global  basis  during  the  month  of 
July  (northern  hemisphere  summer)  than  during  any  other  month  of  the  year. 


Contour  map  of  the  global  representation  of  the  median  value  of  foF2  derived  from 
the  new  coefficients  for  January  1976  at  0600  hrs  UT. 


he  global  representation  of  the  median  value  of  foF2 
ents  for  January  1976  at  1800  hrs  UT. 


Contour  map  of  the  global  representation  of  the  median  value  of  foF2  derived  from 
the  new  coefficients  for  April  1976  at  0600  hrs  UT. 


L» 


Contour  map  of  the  global  representation  of  the  median  value  of  foF2 
the  new  coefficients  for  April  1976  at  1800  hrs  UT. 


Contour  map  of  the  global  representation  of  the  median  value  of  foF2  derived  from 
the  new  coefficients  for  July  1975  at  0600  hrs  UT. 


Contour  map  of  the  global  representation  of  the  median  value  of  foF2  derived  from 
the  new  coefficients  for  October  1975  at  0600  hrs  UT. 


Contour  map  of  the  global  representation  of  the  median  value  of  foF2  derived  from 
the  new  coefficients  for  October  1975  at  1800  hrs  UT. 


tour  map  of  the  global  representation  of  the  median  value  of  foF2  derived  from 
new  coefficients  for  January  1979  at  0600  hrs  UT. 


Contour  map  of  the  global  representation  of  the  median  value  of  foF2  derived  from 
the  new  coefficients  for  January  1979  at  1800  hrs  UT. 


Contour  map  of  the  global  representation  of  the  median  value  of  foF2  derived  fr 
the  new  coefficients  for  April  1979  at  0600  hrs  UT. 


Contour  map  of  the  global  representation  of  the  median  value  of  foF2  derived  from 
the  new  coefficients  for  April  1979  at  1800  hrs  UT. 


3  3 


Contour  map  of  the  global  representation  of  the  median  value  of  foF2  derived  from 
the  new  coefficients  for  July  1978  at  0600  hrs  UT. 


Contour  map  of  the  global  representation  of  the  median  value  of  foF2  derived  from 
the  new  coefficients  for  July  1978  at  1800  hrs  UT. 


Contour  map  of  the  global  representation  of  the  median  value  of  foF2  derived  from 
the  new  coefficients  for  October  1978  at  0600  hrs  UT. 


Contour  map  of  the  global  representation  of  the  median  value  of  foF2  derived  from 
the  new  coefficients  for  October  1978  at  1800  hrs  UT* 
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